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ABSTRACT 

Masses of clusters of galaxies front weak gravitational lensing analyses of ever larger 
samples are increasingly used as the reference to which baryonic scaling relations are 
compared. In this paper we revisit the analysis of a sample of 50 clusters studied as 
part of the Canadian Cluster Comparison Project. We examine the key sources of sys¬ 
tematic error in cluster masses. We quantify the robustness of our shape measurements 
and calibrate our algorithm empirically using extensive image simulations. The source 
redshift distribution is revised using the latest state-of-the-art photometric redshift 
catalogs that include new deep near-infrared observations. Nonetheless we find that 
the uncertainty in the determination of photometric redshifts is the largest source of 
systematic error for our mass estimates. We use our updated masses to determine 
b , the bias in the hydrostatic mass, for the clusters detected by Planck. Our results 
suggest 1 — b = 0.76 ± 0.05(stat) ± 0.06(syst), which does not resolve the tension with 
the measurements from the primary cosmic microwave background. 

Key words: cosmology: observations — dark matter — gravitational lensing — galax¬ 
ies: clusters 


1 INTRODUCTION 


The observed number density of clusters of galaxies as a 
function of mass and redshift depends sensitively on the ex¬ 
pansion history of the Universe and the initial conditions 
of the density fluctuations. Comparison with predictions 
from a model of structure formation can thus constrain cos¬ 
mological parameters, such as the mean density and 
the normalizat i on of the ma t ter p o wer spectrum as (e.g. 
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Henrvl l200dl : iReiprich fc Bohringefl 


20021 ; iHenrv et alilioog) . o r the dark energy equa t ion-o f- 
state w (e.g. Vikhlinin et i _al. 20091: Mantz et al.ll2010l . 120151) . 


For a recent review see lAllen et al.l (120111 ). 


The fact that the observations do not provide actual 
cluster counts as a function of mass, but rather the num¬ 
ber density of objects with certain observational properties, 
such as the number of red galaxies or the X-ray flux within 
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a given aperture, complicates a direct comparison with pre¬ 
dictions: the cosmological interpretation requires knowledge 
of the selection function and the scaling relation between 
the observable and the underlying mass. Furthermore, scal¬ 
ing relations typically have intrinsic scatter that also needs 
to be measured, or at least accounted for. 

One way forward is to simulate the observable proper¬ 
ties of clusters, but the complex non-linear physics involved 
limits the fidelity of such approaches, at least for the mo¬ 
ment. Therefore direct estimates of the clusters masses are 
needed. This can be achieved through dynamical analyses, 
such as the measurement of the motion of cluster members, 
or by measuring the temperature of the hot intra-cluster 
medium (ICM). However, in both cases the cluster is as¬ 
sumed to be in equilibrium, which is generally not a valid 
assumption. For instance, simulations s uggest that hydro¬ 
static X-ray mas s es are biased lo w (e.g. iRasia et al.l |2()0~g1 ; 
iNagai et al'll2007l : lLau et akll2009l ) . 

A more direct probe of the (dark) matter distribution 
would be preferable, which is provided by the gravitational 
lensing distortion of background galaxies: the gravitational 
potential of the cluster perturbs the paths of photons emit¬ 
ted by these distant galaxies, resulting in a slight, but mea- 
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surable, coherent distortion. This in turn provides a direct 
measurement of the gravitational tidal field, which can be 
used to directly infer the projected mass distribution. Note, 
however, that the comparison to baryonic tracers does typ¬ 
ically depend on the assumed geometry of the cluster. For 
a recent review of the use of gravitational lensing to study 
cluster masses and d ensity profiles we refer the reader to 
iHoekstra et al.l (120131 ). 

The sizes of cluster samples are increasing rapidly 
thanks to wide-angle surveys at various wavelengths, es¬ 
pecially at millimeter wavelengths thanks to the inverse- 
Compton scattering of cosmic microwave background 
(CMB) photons off hot electrons in th e intracluster medium, 
the Sunyaev-Zel’dovich Effect (SZE; ISunvaev fe Zeldovichl 
1 1972| j. In fact, the lack of calibrated scaling relations is cur¬ 
rently the dominant limitation of blind searches that exploit 
this effect, such as those carried out using the South Pole 
Telescope (SPT: lReichardt e t all2013il or the Atacama Cos¬ 
mology Telescope (ACT: iHasselfield et ail 1201 -'ll) . The im¬ 
portance of accurate mass calibration is furthermore high¬ 
lighted by the tension between the cosmological parame- 
ters determined from the primary C MB measured by Planck 
llPlanck Collaboration_jrt_alj 2014b|j and those inferred from 


the cluster counts ( Planck Collaboration et al.ll2014c] ). 

Fortunately it is not necessary to determine masses 
for all clusters, but instead it is sufficient to calibrate 
the appropriate scaling relation and its scatter. However, 
doing so still requires substantial samples of clusters for 
which weak lensing masses need to be determined. Even 
for the most massive clusters the uncertainty in the pro¬ 
jected mass is ~ 10%. The triaxial nature of cluster ha¬ 
los, however, leads to an additional intrinsic scatter of ~ 
15 — 20% (e.g.lCorless fc Kindl2007l : iMeneehetti et al.ll20ld : 

iBecker fe KravtsovIhOllI ). Hence to calibrate the normaliza¬ 

tion of a scaling relation to a few percent requires a sample 
of 50 or more clusters. 

To examine the relation between the baryonic proper¬ 
ties of clusters and the underlying matter distribution the 
Canadian Cluster Comparison Project (CCCP) started with 
the study of ar chival observatio ns o f 20 clusters of galax - 
ies, described in lHoekstral ll2007t) and lMahdavi et ahl d2008f> . 
This sample was augmented by observations of an additional 
30 clusters with 0.15 < z < 0.55 with the Canada-France- 
Hawaii Telescope (CFH T). A d etailed descri ption of the 
sample ca n be found in IHoekstra et alj ()2012h (H12 here¬ 
after) and iMahdavi et alJ ( 20131). The compar i son t o the 
X-ray properti e s, pre sented in iMahdavi et alJ (hOOSl) and 
IMahdavi et al.l fl2013h confirmed the predictions from nu¬ 
merical simulations that the hydrostatic mass estimates are 
biased low. 

Other groups have carried out similar studies. The 
Local Cluster Substructure Survey (LoCuSS) used the 
Subaru telescope to carry out a weak lensing study of 

50 clusters, with the most recent results presented in 
lOkabe et alj (120131 ). A thorough analysis of a sample of 

51 clusters was presented by Weighing the Giants (WtG; 
Ivon der Linden et al.l l2014al : lApplegate et al.l I20l4) . For a 
large fraction of the clusters the latter s tudy also obtaine d 
photometric r edshifts for the sou rces d Kelly et al.l 12014) . 
Most recently lUmetsu et al.l d20l4 ) presented results for a 
sample of 20 massive clusters, 17 of which were observed by 
WtG. In general there is significant overlap as these stud¬ 


ies all target ma ssive well-known cluste rs of galaxies. This 
was exploited bv lApplegate et~aTI (12014 ) who compared the 
masses from the various studies. Although they find an ex¬ 
cellent correlation with the results from H12, the CCCP 
masses are on average ~ 20% lower than their estimates. 
This is much larger than the statistical uncertainties and 
warrants further investigation. This is the aim of this paper. 

A correct interpretation of the inferred weak lensing sig¬ 
nal relies on accurate shape measurements and knowledge of 
the redshifts of the sources used in the analysis. The former 
has been examined quite extensively over the past decade, 
for instance in severa l blind studies using simulated images 
dHevmans et al.ll2006l:lMassev et aTlhOOTlJBridle et al.lhOld : 

iKitching et al.ll2012l : lMandelbaum et al. l2014a| j . The results 

of such simulations have been used to quantify the biases in 
shape measurements, but the sensitivity of the calibration 
to the input of the simulations has not been investigated in 
much detail. However, thanks to an improved understand- 
ing of the sources of bias, and how they propagate (e.g . 

Massey et al.|20Lj|: ISemboloni et al.ll2013l : iMiller et al.l2013 : 

Viola et al.lhOldfh it has become evident that a correct inter¬ 

pretation of these simulations depends critically on how well 
they match the specihc observations under consideration. In 
§2 we examine the importance of the fidelity of the image 
simulations. To calibrate our method, we create an extensive 
set of images, varying a number of input parameters. 

Another important source of uncertainty is the redshift 
distribution of the sources. In §3 we present our photomet¬ 
ric redshift estimate s based on measurements in 29 bands in 
the COSMOS held (IScoville et al.l [20071 : ICapak et al1l2007f) 
including n ew deep observations in five NIR bands from Ul- 
traVISTA (iMcCracken et akll2012h . We also revisit the issue 
of contamination by cluster members. We present new weak 
lensing mass estimates in §4 an d use these in §5 to calibrate 
t he h ydrostatic masses used bv IPlanck Collaboration et al.l 
(l2014ch to infer cosmological parameters. Throughout the 
paper we assume a cosmology with = 0.3, I2 a = 0.7 and 
Hq = 70/i7o km/s/Mpc. 


2 CALIBRATION OF SHAPE 
MEASUREMENTS 

The measurement of the shapes of small, faint galaxies is 
one of two critical steps in order to derive accurate clus¬ 
ter masses from weak gravitational lensing, the other step 
involving knowledge of the source redshift distribution. We 
discuss the latter in S|3] and focus first on the algorithms 
used to measure galaxy shapes. Most studies to date have 
focussed on the correction for the blurring by the PSF, which 
leads to rounder images (due to the size of the PSF) and pre¬ 
ferred orientations (if the PSF is anisotropic). An incomplete 
correction for the former leads to a multiplicative bias fj, and 
a residual in the latter to an additive bias c; the observed 
shear and true shear are thus related by (e.g. lHevmans et al.l 

hood) : 

7 ? bs = (l + /i) 7 r ue + c, (1) 

where we implicitly assumed that the biases are the same 
for both shear components. For cosmic shear studies the ad¬ 
ditive bias is a major source of concern because t he (resid - 
ual) PSF introduces power on relevant scales (e.g. iHoekstral 
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120041 ). For cluster lensing the additive bias is less impor¬ 
tant because the measurement of cluster masses involves 
the azimutlially averaged tangential shear and PSF patterns 
largely average out for our data. We study the residual addi¬ 
tive bias in Appendix El and find that we can indeed ignore 
the residuals arising from PSF anisotropy in our analysis. 

One approach to recover the true galaxy shape is to 
assume a suitable model for the galaxy light distribution, 
which is subsequently sheared, convolved with the PSF and 
pixellated. The model parameters are varied until a best fit 
to the data is obtained. This has the advantage that the 
detrimental effects of the PSF (and other instrumental bi- 
ases) can be incorporated into a Bayesian fram ework (e.g. 
iMiller et a.l Jl20 l73 ; iBernsI ein fe Armstronell2014l) . The chal¬ 
lenge, however, is to use a model that provides a good de¬ 
scription of the galaxies, while having a limited number of 
parameters in order to avoid over-f itting. A model that is too 
rigid will lead to model bias (e.g. iBernsteinl 12010 ). whereas 
a m odel that is too flexib le tends to fit noise in the images 
fe.g. lKacprzak et al.ll20l3) . Furthermore, accurate priors for 
the size and cllipticity distributions (and any other param¬ 
eter entering the model) are required to obtain an unbiased 
estimate for the shear. 

An alternative approach, which we use here, involves 
measuring the moments of the galaxy images, which are sub¬ 
sequently corrected for the PSF. The shapes can be quanti¬ 
fied by the polarization: 


n — 12 2 . 2/12 

In + I 22 ’ 62 “ /11 + / 22 ’ 


( 2 ) 


where the quadrupole moments Iij are given by 


1 

To 


d 2 xii 2 ; j W (x) / (x), 


(3) 


where /(x) is the observed galaxy image, VF(x) a suitable 
weight function to suppress the noise and Io the weighted 
monopole moment. I 11 the case of unweighted moments, Io 
corresponds to the flux and the correction for the PSF is 
straightforward as the PSF corrected moments are given b£l 

7-true r-obs xPSF (A\ 

-*-ij ■‘■ij ? VV 


i.e., one only needs to subtract the moments of the PSF from 
the observed moments. The result provides an unbiased es¬ 
timate of the polarization. However, the change in polariza¬ 
tion 5a due to a shear depends on the unsheared shape 
eint: it is more difficult to change the shape of an object 
that is already elongated. This response is quantified by the 
polarizability P 7 , defined such that 5d = P^S'yi. As the 
shear is obtained from an ensemble of galaxies, an unbiased 
estimate thus req uires knowledge of the intrinsic ellipticity 
distribution ('e.g. lviola et al.ll2014l ). 

Unfortunately real data contain noise and thus un¬ 
weighted moments are not practical. To suppress the ef¬ 
fects of noise a weight function needs to be chosen, ide¬ 
ally matched to the size and s hape of the galaxy im- 
age. However, as discu ssed in e.g. iMassev et al.l (2013) and 
ISemboloni et alj (120131 ) this complicates matters as the cor¬ 
rection for the PSF now involves higher order moments, 


1 We assume that the measurement is centered on the location 
where the dipoles vanish. 


which themselves are affected by noise. Limiting the expan¬ 
sion in moments is similar to the model bias in fitting meth¬ 
ods. 


Recent studies using simulated data have shown that 
multiplicative bia ses depend strongly on the signal-to- 
noise ratio (SNR ; iBridle et al.l l2010l : iKitching et al.l 120121 : 
IMiller et al.ll20li)l ) with some hints alre ady present in the de - 
pendence of the bias on magnitude in IMassev et al.l d2007t) . 
As the origin of this bias is now better understood, it has 
also become clear that the performance of a particular al¬ 
gorithm will depend on the data it is applied to. Hence the 
performance evaluation, such as the determination of the 
bias that one wishes to correct for, depends on the input of 
the simulations: if the input does not match the actual data, 
the inferred bias may be different from the actual value. Al¬ 
though SNR is the most critical parameter, the bias may 
also depend on th e galaxy profile, or the size and ellipticity 
distributions (e. g. Miller et al .1 20131 : iMelchior fe Violall2012l ; 
iKacprzak et al.ll2012l : Viola et al.l 20141) . Unless the fidelity 
of the simulation can be somehow guaranteed, the sensitivity 
of a method to the input parameters needs to be quantified 
and the uncertainties propagated. 

In this paper we focus on the commonly used 
KSB method de v eloped by iKaiser et alJ (1 19951 ) and 
Luppino fe Kaiser! (1 19971 ) with correction s pro vided in 


Hoekstra et al. ( 199sl) and iHoekstra et al.l (12000). It was 
used t o determine masses for the CCCP sample in Hoekstral 
(|2007t) and H12; we refer the interested reader to these pa¬ 
pers for more details. The object det ection is done usin g 
the hierarchical peak finder described in lKaiser et al.l(jl995l ). 
which gives an estimate for r g , the Gaussian scale radius of 
the object. This value is used to co mpute the weighted mo - 
ments, which are corrected following IHoekstra et al.l (1 9981 ). 
In addition we also compute a e , the uncertainty in the polar¬ 
ization, which is approximately oc 1 /u, where v is the signal- 
to-noise ratio of the detection (IHoekstra et al.ll2000l ). This 
allows us to downweight the noisy galaxies and we therefore 
estimate the average shear for an ensemble of galaxies as 


hi) 


E W i 


with Wi 




(5) 


where (e 2 ) is the intrinsic variance of the galaxy elliptic¬ 
ity components. This is the dominant source of uncertainty 
for the shear for bright objects, a nd we adopt a value of 
(e 2 ) 1 / 2 = 0.25 l|Hoekstra et al.l200d ). In our image simations 
we vary the input cllipticity distribution (see 3231), which 
in principle would require adjusting the value for (e 2 ) 1 ^ 2 
accordingly to optimally weight objects. However, for sim¬ 
plicity we keep it fixed when we quantify the multiplicative 
bias. 


2.1 Input galaxy properties 

To populate our image simulations we use a sample of galax¬ 
ies for which morphological parameters were meas ured from 
resolv ed F606W images from the GEMS survey (jR.ix et al.l 
120041) . These ga laxies were mode led as single Sersic mod¬ 
els with galf it dPeng et al.ll200j) and for our study we use 
the half-light radius, magnitude and Sersic index n. We only 
consider galaxies fainter than m r = 20 because bright ob¬ 
jects might cause unrealistic features in the simulated im- 
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Figure 1. Number counts as a function of r-band magnitude for 
our image simulations (solid black histogram), STEP1 (dotted 
red histogram) and STEP2 (dashed blue histogram). The counts 
were normalized in the range 20 < m r < 22 and we adjusted 
the magnitudes of the STEP simulations for the use of different 
filters. 


Figure 2. Input ellipticity distributions, described by a Rayleigh 
distribution truncated at e = 0.9, for eo = 0.35 (thick black 
curve), eo = 0-1 (reel dotted curve) and eo = 0.5 (blue dashed 
curve). For reference we also show a histogram of ellipticities from 
one of the STEP1 simulations. 


ages. Excluding these does not impact our results as we do 
not use them in our source sample anyway: to measure the 
lensing signal we use galaxies with 22 < m r < 25. 

The resulting number density of galaxies as a function 
of apparent magnitude m r is presented in Figure |T] (solid 
black histogram). The results suggest a power-law for the 
counts, where the flattening for m r > 25.5 is caused by 
incompleteness of the input catalog. In principle, faint un¬ 
resolved galaxies can affect shape measurements of brighter 
galaxies through modulation of the background noise and 
blending. In 1 12.31 we therefore examine the need to include 
fainter galaxies in the simulations. 

H12 based their assessment of the accuracy of the shape 
measurements on the results from the Shear TEsting Pro - 
grammes (STEP: [Hevmans et al.ll2006l ; iMassev et ai1l2007l) . 
These were blind challenges with the aim to benchmark the 
performance of shape measurement algorithms, especially 
for cosmic shear studies. In both cases the implementation 
used by H12 performed well, with an average multiplicative 
bias of ~ 2%. As a consequence, H12 ignored the multiplica¬ 
tive bias in t heir mass estimates. 

STEP1 dHevmans et al.ll2006l ) simulated CFHT obser¬ 
vations in the /-band with an integration time of 3600s, 
which should be quite comparable to our data (a total in¬ 
tegration time of 1 hour in the r' band using CFHT). The 
red dotted histogram in Figure [T] shows the galaxy num¬ 
ber counts that were used as input for STEP1 (converted 
to r-band assuming a mean galaxy color of r — i = 1). The 
counts are normalized such that the sum is the same for all 
three examples in the range 20 < m r < 22. Comparison to 
the GEMS catalog shows that STEP1 lacks the faint galax¬ 
ies that are present in real observations, even if they are 


unresolved. As we will show in £12.31 this leads to a signifi¬ 
cant underestimate of the multiplicative bias for the actual 
CCCP data. _ _ 

STEP2 (IMassev et al.l 120071) simulated images that 
would be obtained with an exposure time of 40 minutes 
in good conditions with SuprimeCam on Subaru. Given 
the larger aperture and throughput of Subaru compared to 
CFHT this corresponds to a total exposure time that is ~ 4 
times longer than the CCCP data. The input galaxy im¬ 
ages were based on a shapelet decomposition of resolved 
galaxy images from the Hu bble Space Telesc ope COSMOS 
survey JScoville et al.l 12007I ). as described in IMassev et al.l 
(12004 ). As a result the simulations should better capture the 
complex morphologies of real galaxies. The number counts, 
shown by the blue dashed histogram in Figure [l] match the 
GEMS input counts much better than STEP1, although in¬ 
completeness occurs at m r ~ 24.5. 

In addition to the magnitudes, the GEMS catalog pro¬ 
vides values for the Sersic index of the galaxies, as well as 
their half-light radii and ellipticities. The use of Sersic pro¬ 
files to descri be the g a laxies may limit the fidelity of the sim¬ 
ulations (see iKacprzak et al.ll2014l . for a study of the biases 
that may arise). We examine the bias as a function Sersic 
index in USD and find that our shape measurement algo¬ 
rithm is not particularly sensitive to the profile, especially 
when compared to other sources of bias. We therefore expect 
that the difference with using realistic galaxy morphologies 
is small. One of the aims of the third g ravitationa l lensi ng 
accura cy testing challeng43 (GREAT3; iMandelbaum et alJ 
l2014bl) is to compare the results of shape measurement 


2 http://www.great3challenge.info 
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methods on postage stamps of actual (PSF corrected) HST 
observations and the corresponding Sersic fits. 

The different parameters describing the galaxies are 
jointly sampled from the GEMS catalog to account for their 
intrinsic correlations (e.g., brighter galaxies are on average 
lar ger). How e ver, w e do not use the ellipticities provided 
by iRix et al.l (120041') . because of concerns that these do not 
match our data, as discussed in £12.21 Instead we use a para¬ 
metric description, which allows us to investigate the role of 
the ellipticity distribution. We assign ellipticit ie43 6 that are 
drawn from a Rayleigh distribution given by 


-P(e,eo) 



(6) 


where the value of eo determines the width of the distribu¬ 
tion, as well as the average (e) = eo s/tT/2. We need to trun¬ 
cate the distribution because the ellipticity cannot exceed 
unity, but also because galaxy disks have a finite thickness. 
We therefore set P(e, eo) = 0 if e > 0.9. We assume that the 
ellipticity distributi on is independen t of ot her galaxy prop¬ 
erties, whereas e.g.. Ivan Uitert et alj (l201l| j did observe dif¬ 
ferent distributions for early and late type galaxies. Given 
the accuracy we require here, we find that this assumption 
does not impact our results. The ellipticity distribution of 
the GEMS catalog matches that of eo = 0.35 for moderate 
ellipticities dMelchior fc Violal 120121 ') which is indicated by 
the black curve in Figure 0 We also show input ellipticity 
distributions for eo = 0.1 (red dotted curve), and eo = 0.5 
(blue dashed curve). For comparison we al so show the in¬ 
put e llipticity distribution used by STEP1 dHevmans et alj 
120061 ). which peaks at very low ellipticities. 
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Figure 3. Comparison of the simulations and the actual CCCP 
data for Abell 1835 (black lines) for bright and faint sources. 
The left panels show the observed polarizations, i.e. uncorrected 
for the PSF. We find that the ellipticity distribution for bright 
galaxies (20 < m r < 22) is best approximated if we take eo = 
0.15, whereas the ellipticity distribution for fainter galaxies (22 < 
m r < 25) requires a higher value of eo = 0.25. For reference we 
show eo = 0.35 which is ruled out by the data. As shown in the 
right panels, the simulations match the observed distribution of 
half-light radii fairly well. 


2.2 Description of the simulations 

To create the images we use galsim (|Rowe et al.l 120141 ), a 
publicly available code tha t was developed for GREAT3 
iMandelbaum et al.l l2014bh . The main input is a list of 
galaxies with a position, flux, lialf-light radius, Sercic index 
and ellipticity, from which sheared images are computed. We 
limit the sample to objects with 0.5 < n < 4.2 because of 
limitations of the version of galsim we used. To create the 
simulated images we draw objects from the GEMS catalog. 
Given the limited number of galaxies observed by GEMS, 
objects typically appear multiple times in the simulation, 
but with a different ellipticity and orientation. 

The weak lensing analysis of CCCP Megacam data is 
done on stacks with a total integration time of 1 hour each. 
For each cluster two such stacks are observed, which are 
merged at the catalog stage (see H12 for details). Our im¬ 
age simulations therefore assume the same noise level as ob¬ 
served in these data. To simulate the observed data we also 
need to provide a realistic PSF, for which we use a circu¬ 
lar Moffat profile with a FWHM=0 / . , 67 and /3 = 3.5. This 
resembles our observations of Abell 1835, which are typical 
for the CCCP sample, and we also adopt the noise level ob¬ 
served in these data. We include a low number of stars to 
measure the PSF in the images. In Appendix [B] we quantify 

3 The ellipticity is defined as (a — b )/ (a + b ), with a and b the 
major and minor axes, resp. The polarization for such a galaxy 
would be ~ (a 2 — b 2 )/(a 2 + b 2 ) 


the impact of realistic star densities. We find that the ob¬ 
served star densities in the CCCP data are sufficiently low 
that they do not impact the results. 

We create pairs of images where the galaxies are rotated 
by 90 degrees in the second image to reduce the noise due 
to th e intrinsic ellipticity distribution (see e.g. iMassev et al.l 
l2005t) : by construction the mean intrinsic ellipticity when 
both are combined is zero. We analyse the images separately 
and thus, due to noise in the images, this is no longer exactly 
true, especially for faint galaxies. The input shears typically 
range from -0.06 to 0.06 in steps of 0.01 (for both compo¬ 
nents), yielding 169 image pairs for each ellipticity distribu¬ 
tion. Each image has a size of 10,000 by 10,000 pixels, with 
a pixel scale of 0V185, the same as our MegaCam data. This 
results in a sample of ~ 10 7 galaxies with 20 < m r < 25 for 
each value of eo- To examine the dependence on seeing and 
PSF anisotropy we create somewhat smaller sets, consisting 
of 49 pairs of images. 

We analyse these images in the same way as the CCCP 
data. Figure [3] shows the distribution of observed polar¬ 
izations and half-light radii from the actual data (solid 
histograms) and simulated data (dotted and dashed his¬ 
tograms). We reproduce the magnitude distribution (not 
shown) and the size distribution for galaxies fainter than 
m r = 22. As shown in the top right panel of Figure [3] there 
are many simulated bright galaxies that have half-light radii 
that are large compared to our CCCP data. For these galax¬ 
ies the polarizations are significantly smaller than the dis- 
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Figure 4. Multiplicative bias for sources with 20 < m r < 25 as a 
function of eo (see text). The red squares show the results for the 
implementation of KSB used in H12. The black points show the 
results when we use the observed estimate for P 7 for individual 
galaxies. Although the bias is larger in the latter case, it depends 
only weakly on eo, especially in the relevant range of 0.15-0.3. 
The open circles indicate the bias if SExtractor is used instead 
of hf indpeaks. 

tribution with eo = 0.35 found bv lMelchior fc Violal (|2012l ). 
Although the input catalog is based on HST data, we suspect 
that the use of galf it may give too much weight to the outer 
regions of the galaxies, which are down-weighted in moment- 
based methods. This highlights the difficulty in establishing 
the input ellipticity distribution, which remains rather un¬ 
certain. For the main sample of sources, with 22 < m r < 25, 
we find a good match if we adopt eo = 0.25. We take this 
value as our reference in the remainder of this paper. We 
will conservatively assume that 0.15 < eo < 0.3 when we 
estimate systematic uncertainties in the empirical bias cor¬ 
rection. 

2.3 Multiplicative bias as a function input 
parameters 

As the underlying ellipticity distribution remains uncertain, 
we start by examining the average bias of KSB as a func¬ 
tion of eo for galaxies with 20 < m r < 25, i.e. the range 
in magnitude of the sources used in the CCCP analysis by 
H12. To det ect objects we use t he hierarchical peak finder 
described in iKaiser et al.l (1 19951 ). which is the default al¬ 
gorithm in our analysis. The main difference between the 
various implementations of the KSB algorithm is the way 
the shear polarizability P 1 is estimated. As the observed 
values are noisy, H12 used a parametric fit to average val¬ 
ues as a function of size for different magnitude bins (also 
see lHevmans et aI1l2006l . for a concise description). The red 
squares in Figure[4]show the bias as a function of eo for this 
implementation of KSB. The bias changes by 0.04, which 


Figure 5. Multiplicative bias as a function of apparent magni¬ 
tude for the simulated CCCP data. The measurements are av¬ 
erages for bins with a width of 0.5 magnitude. The black points 
correspond to an elliptictity distribution with eo = 0.25. The red 
and blue lines are for eo = 0.15 and eo = 0.3, respectively. The 
red squares indicate the bias if we follow the procedure to evalu¬ 
ate P 7 used in H12. For comparison, the hatched region indicates 
the 68% confidence interval for the average bias for galaxies with 
20 < m r < 25 (for e 0 = 0.25). 

corresponds to a relative change of ~ 40%, over the rather 
extreme range in ellipticity distribution. For eo = 0.25 we 
find a bias of fj, ~ —0.115, which is mu ch la rger than the 
value reported in iHevmans et al.l (120061 ) and iMassev et al.l 
(|200 7t ). Consequently we cannot ignore the multiplicative 
bias, as was done in H12. 

The black points show the results if we use the mea¬ 
sured value of P 7 for each galaxy. In this case the bias is 
larger (fi ~ —0.165 for to = 0.25), but also less sensitive 
to the ellipticity distribution. For this reason, as well as 
simplicity, we adopt this implementation as our reference. 
We note that the full chain of detection and shape analy¬ 
sis needs to be simulated. This is highlighted by the open 
black points, which indi cate the bias if we use SExtractor 
dBertin fe Arnoutsll 19961 ) to detect objects and use the value 
for FLUX_RADIUS to compute the corresponding value for r g : 
the observed bias is affected at the per cent leveJJ 

Figure [5] shows that the bias increases quickly for fainter 
galaxies, irrespective of the ellipticity distribution. This is 
also the case when we consider the implementation used 
by H12 (red squares). A strong dep endence of th e bias 
on the SNR was already observed in iBridle et al.l (l20ld ) 

4 In the process of making this comparison we discovered that 
SExtractor (we used version 2.5.0) incorrectly assigns objects 
FLAG=16 if they are elongated horizontally. This problem can be 
avoided by adopting a value for MEMORY-BUFSIZE larg er than the 
image dimensions. We note that iGruen et alj ||2014| ) discovered 
the same problem and reported on this. 
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Figure 6. Multiplicative bias as a function of seeing for galaxies 
with 20 < m r < 25 for an elliptictity distribution with eo = 0 
(red line), the reference case with eo = 0-25 (black points) and an 
extreme case with eo = 0.5 (blue line). The histogram shows the 
distribution of PSF sizes of the CCCP data measured for each 
chip. The image quality is typically best in the inner regions of 
the field-of-view, which are most relevant for the mass estimates. 

and iKitching et al.l (|2012l ). The lack of faint galaxies in 
S TEP1 is the main reas on that a small bias was observed 
in lHevmans et all (|2006l ). When we restrict the analysis to 
the magnitude range simulated by STEP1 we reproduce the 
small bias for the implementation used for that paper. Note 
that STEP2 simulated data that are deeper than our CCCP 
data. The implementation used by H12 gives smaller biases 
when considering the full range in magnitude, but overcor¬ 
rects bright galaxies (i.e. /r > 0). It appears that the choice 
of the fitting function partly compensated for the bias due 
to noise. 

The SNR. is also affected by the PSF size: the larger 
the PSF, the lower the SNR as the flux is spread over more 
pixels. The seeing also determines how well galaxies are re¬ 
solved, which impacts the bias as well (see Appendix 0. 
Figure[5]shows the value of fi for galaxies with 20 < m r < 25 
as a function of seeing for eo = 0 (red line), eo = 0.25 (black 
points) and an extreme case with eo = 0.5 (blue line). Note 
that we keep the range in apparent magnitude the same. The 
results demonstrate the importance of good image quality: 
the bias more than doubles from -0.11 to -0.25 as the seeing 
deteriorates from O'/5 to 1". 

The number of faint galaxies increases rapidly (cf. 
Fig. CD, which results in source galaxies blending with 
fainter ones. Even if a faint galaxy is not detectable, it 
will impact the noise level, effectively introducing corre¬ 
lated noise that affects the local background determina¬ 
tion. Both of these will modify the multiplicative bias in 
a way tha t that can o nly b e quantified through simu¬ 
lations. In iBridle et al.l (I2010I ). IKitching et al.l (| 20121 ) and 
iMandelbaum et al.l ( 2014ah only postage stamps of isolated 


Figure 7. Multiplicative bias for an ellipticity distribution with 
eo = 0.25 where only galaxies with magnitudes brighter than mn m 
are included in the simulation. The black points show the bias for 
20 < m r < 25. For comparison, the hatched region indicates the 
68% confidence interval for the average bias for galaxies using the 
GEMS input catalog. Irrespective of the magnitude range, the 
bias converges when mn m is 1.5 mag fainter than the magnitude 
limit of the source sample. 


galaxies were analyzed, and thus the effects of blending and 
faint galaxies were not included. Figure [3 shows that this 
is an important effect, and cannot be neglected. To obtain 
these results we create images where we include galaxies 
down to a limiting magnitude mii m . The input GEMS cata¬ 
log is incomplete for m r > 25.5 (see Figure Q} and we aug¬ 
ment the catalog by duplicating the fainter galaxies such 
that the input counts follow the power-law relation seen at 
brighter magnitudes. At the faintest magnitudes these galax¬ 
ies are unresolved in the simulated ground-based data, and 
hence the details of their structural properties are not criti¬ 
cal. 

The black points in Figure [3 show that the bias for 
galaxies with 20 < m r < 25 increases until mn m > 26.5; in 
general we find that the bias converges if we include sources 
that are 1.5 magnitude fainter than the magnitude limit 
of the sample of sources used to measure the weak lensing 
signal. This also appears to be true if we consider narrow 
bins in magnitude, such as the bin with 24 < m r < 25 for 
which the bias is large, but converges for m r > 26.5. The 
dominant contribution of these faint galaxies is to act as a 
source of correlated noise, affecting the shape measurements 
of brighter galaxies. These results demonstrate that it is 
important to ensure that the input catalog used for image 
simulations contains a sufficient number of galaxies fainter 
than the magnitude limit one is interested in. 

For comparison the hatched area in Figure [3 indicates 
the 68% confidence region for the bias we obtain when we 
use the GEMS input catalog, without introducing additional 
faint galaxies to account for incompleteness (for m r > 25.5). 
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Figure 8. Multiplicative bias as a function of input ellipticity dis¬ 
tribution for four different Sersic indices (as indicated) for galaxies 
with 20 < m r < 25. The hatched region indicates the 68% confi¬ 
dence region when the distribution of Sersic indices from GEMS 
is used. The bias depends on the value of the Sersic index, al¬ 
though we note that the results cannot be compared directly, as 
explained in the text. 


Comparison with the black points suggests that the input 
catalog is sufficient for the interpretation of the CCCP data, 
and we use it to compute the corrections in 1)2.41 

Finally we examine whether the bias depends on the 
assumed distribution of Sersic indices. To do so, we create 
images where all galaxies have the same Sersic index n, while 
keeping the other parameters the same. The results are pre¬ 
sented in Figure[8]for different values of eo- The bias depends 
on the value of n, although the range is small (~ 0.02 for 
eo = 0.3), with the results for n = 1 (corresponding to expo¬ 
nential profiles) most discrepant. We note, however, that the 
half-light radii were kept to the values listed in the GEMS 
catalog, which can lead to small changes in the SNR, com¬ 
plicating a direct comparison. Given the small variation in 
/x, and the fact that the observed distribution of Sersic in¬ 
dices is well constrained, we assume that the uncertainty in 
this distribution can be ignored. Hence, the dominant un¬ 
certainty in our bias estimate arises from the uncertainty in 
the cllipticty distribution. 


2.4 Empirical correction 

If the simulated data resemble the actual observations suffi¬ 
ciently well, the average bias for the source sample could be 
used to adjust the cluster masses accordingly. Note that one 
would still have to determine the bias as a function of see¬ 
ing. An additional complication arises, however, because we 
lack redshifts for the sources: the bias depends on the fluxes 
and sizes of the sources. As a result the bias may be redshift 
dependent, which is not captured by the image simulations 
as the same shear is applied to all sources. In reality the am¬ 


plitude of the shear signal depends on the geometry of the 
lens-source config uration, quantified by the criti cal surface 
density E cr i t (e.g. iBartelmann fc Schneideiil200ll ) 

v . _ c D* f7\ 

crlt 4 t vGDtDu’ U 

where D s , Di and Di s are the angular diameter distances 
between the observer and the source, the observer and the 
lens, and the lens and source, respectively. The sensitivity 
to the source redshift distribution is quantified by the ratio 
P = Di s /D s . The average shear for an ensemble of galaxies 
is proportional to ((1 + /x)/3). If photometric redshifts for 
the individual sources are available, the redshift dependence 
of the lensing signal can be accounted for on an object-by- 
object basis and an average correction for the multiplicative 
bias is possible. An alternative route, which we take here, is 
to compute the multiplicative bias using the observed prop¬ 
erties of individual galaxies. The correction, however, will 
still depend on the intrinsic ellipticity distribution of the 
sources. 

We assume th at the bias is only a function of signal- 
to-no ise ratio (e.g. I Melchior fc Violal 120121 : iKacnrzak et al.l 
2012i) and the size relative to that of the PSF (e.g. 
Massey et al 120131) . We quantify the latter by the parameter 
7Z , defined as: 


n 2 = 


r h,* 


h, gal 


(8) 


where rh ,* denotes the half-light radius of the PSF and r^.gal 
that of the observed galaxy. Despite being a simple prescrip¬ 
tion, we show in 1 )2.51 that this captures the dependence on 
PSF size quite well. As a proxy for the signal-to-noise ratio 
we take v = 1 /q e , the reciprocal of the uncertainty in the 
polarization dHoekstra et al.M2000l ). We refer the interested 
reader to Appendix □ for more details about our empirical 
correction, which is given by 




1 + a{to)TZ 


(9) 


The dependence on the resolution parameter 1Z is described 
by a single parameter a that is a function of eo only. We 
require three free parameters to describe dependence of the 
bias on v. b(v) = &o + bi/y/(v) + & 2 /with fit parame¬ 
ters that vary smoothly with eo. The best fit parameters 
as a function of eo are listed in Table IG1I Although this 
parametrization does not describe the simulated data per¬ 
fectly, and obvious improvements can be suggested, we pre¬ 
fer our choice as it provides a sufficently accurate correc¬ 
tion, with a relatively small number of parameters. Includ¬ 
ing more parameters did not improve the robustness of the 
correction. 

We find that our parametrization of the bias does not 
perform well for galaxies with large observed sizes (rh > 5 
pixels; see e.g. Figure [Gill . As discussed in more detail in (2] 
we find that the recovered lensing signal for these galaxies is 
biased low (see Figure fl4ll . Closer investigation of the simu¬ 
lated data shows that most of these galaxies are intrinsically 
small and faint. In some of the cases the sizes are increased 
by noise in the images, but a large fraction is blended with 
other galaxies. The large increase in galaxy density in clus¬ 
ters of galaxies is expected to exacerbate this problem, which 
is not captured by our simulations (which are representative 
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Figure 9. Residual multiplicative bias as a function of input ellip- 
ticity distribution for galaxies with 22 < m r < 25 and ry, < 5 pix¬ 
els. The black points indicate the results when the correct value 
for eo is used in the correction. The line indicates the residual bias 
if we assume eo = 0.25 in the correction, instead of the correct 
value for the simulated distribution. Adopting a value eo = 0.25 
for the correction results in // C or < 0.015 over the expected range 
in actual eo values (indicated by the hatched region). 

of the field). To minimize the impact of blending, we include 
only galaxies with rh < 5 pixels in the lensing analysis. This 
size cut is applied to the tests presented below, as well as 
our actual measurements. 

2.5 Testing the empirical correction 

To quantify how well the correction works when we apply 
it to the simulated data, we first examine the residual bias 
/r cor as a function of eo- As explained in more detail in H3.ll 
we restrict the source sample to galaxies with 22 < m r < 25 
to allow for a better correction for the contamination by 
cluster members. In addition we apply a size cut, requiring 
that rh < 5 pixels. This is motivated by our image simula¬ 
tions where we found that the correction for large galaxies 
is biased, because they are blended or too faint to have their 
shapes measured reliably. We therefore limit the discussion 
of the performance of the empirical correction to this range 
in apparent magnitude and galaxy size. The results are pre¬ 
sented in Figure [9] which shows that for the range of interest 
for eo (indicated by the hatched region) |/x cor | < 0.005. 

As the intrinsic ellipticity distribution remains uncer¬ 
tain, it is useful to examine the bias that is introduced when 
an incorrect value for eo is used for the empirical correction. 
If we use the parameters corresponding to eo = 0.25 to cor¬ 
rect the measurements from other input distributions we 
find that /r cor is still small, as indicated by the black line in 
Figure^ Our empirical correction is quite robust against the 
uncertainty in the input ellipticity distribution (if we take 
eo = 0.25). As discussed in Appendix [Cl the parametrization 



seeing [”] 


Figure 10. Residual multiplicative bias as a function of seeing 
for sources with 22 < m r < 25. The black points show the results 
for eo = 0.25. The red (blue) hatched regions indicate the 68% 
confidence region for the bias if we use the parameters for eo = 
0.25 to correct the simulations with input distribution with eo = 
0.15 (eo = 0.3) instead. The bottom histogram shows the seeing 
distribution for each chip in CCCP, the side histogram shows the 
corresponding distribution of residual bias, with ( ficor) = —0.001. 


for the size dependence of the bias is not accurate for large 
galaxies, which are typically bright. This is indeed reflected 
in the residual bias as a function of apparent magnitude: we 
observe /r cor ~ 0.02 for m r < 22, with a bias ~ 0 for galaxies 
with m r > 22. 

The empirical correction was determined for a partic¬ 
ular PSF and integration time. Although it is in principle 
possible to create simulated data sets for each set of observ¬ 
ing conditions, a useful correction scheme should be more 
generally applicable. As discussed in Appendix [D] we also 
simulated data from the second Red-sequence Cluster Sur¬ 
vey (RCS2). These data are shallower, but the results pre¬ 
sented in Figure [Dl] indicate that the correction works fairly 
well for these shallower data. This suggests that the model¬ 
ing of the SNR-dependence is adequate. 

More interesting is whether our approach to quantify 
how well galaxies are resolved, i.e. the choice of TZ, can be 
used for a range of seeing values. To this end we correct the 
set of images used to study the seeing dependence of the 
bias (see Fig. [6]). The results for galaxies with 22 < m r < 25 
are presented in Figure 111)1 which shows /i cor as a function 
of the FHWM of the PSF. Even for a FWHM of 1" the 
bias is reduced significantly. Nonetheless the residual bias 
can still be substantial. However, as is shown by the seeing 
histogram, the CCCP data span a relatively narrow range, 
and the mean bias for the full sample is (/r cor ) = —0.001. 
Furthermore, the largest FWHM values occur on chips far 
away from the cluster location. We therefore ignore the see¬ 
ing dependence, as the residual bias is still much smaller 
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than the statistical uncertainties for individual clusters, and 
the ensemble average bias very small. 

We conclude that our empirical correction is adequate 
to determine cluster masses for the CCCP sample. Based 
on the residuals we assign a systematic uncertainty of 2% 
in the cluster masses due to the uncertainty in the input 
cllipticity distribution, the limited exploration of the role of 
morphology, and the variation in image quality. 


3 SOURCE REDSHIFT DISTRIBUTION 

We lack the color information to derive photometric 
redshifts for th e indiv idual sources, as o p posed to e.g. 
lAppleeate et al.l (120141 ') and lUmetsu et al.l (120141 ). Fortu¬ 
nately it is sufficient to know the average source redshift 
distribution, which we discuss in more detail in We 

note, however, that photometric redshifts enable an opti¬ 
mal weighting of the sources. In particular photometric red¬ 
shifts allow for a separation of source galaxies and cluster 
members. The latter are unlensed, as are foreground galax¬ 
ies, and thus dillute the observed lensing signal by a factor 
/contain (r) = 1 + n c i (r)/rifld, where n c i(r) is the number den¬ 
sity of cluster galaxies and rifld the number density of field 
galaxies. This correction is especially important at small dis¬ 
tances from the cluster center. As we describe in 43.1l we can 
correct for the reduction in signal by quantifying the level 
of contamination. This assumes that the orientations of the 
cluster members are rand om, which is supported by obser¬ 
vations (ISifon et al.ll20lJ l. 

In doing so, we assume that the change in counts is 
solely caused by contamination by cluster members. How¬ 
ever, gravitational lensing not only changes the galaxy 
shapes, but also magnifies sources. As a result the back¬ 
ground sources appear brighter, leading to an increase in 
the observed counts. On the other hand, the actual volume 
is reduced, because the observed solid angle corresponds 
to a smaller solid angle behind the cluster. Consequently, 
the net change depends on the number density of back¬ 
ground gala xies as a fun ction of apparent magnitude (see 
e.g. §3.4 in iMellierl Il999h . In our case we observe a slope 
dlogAgai/dM ~ 0.38 — 0.4 for galaxies with 22 < m r < 24. 
T his is somewhat ste eper than the slope of ~ 0.33 observed 
bv lHogg et al.l (Il997l ) in the R— band. In either case the net 
effect is minimal: even for k = 0.1 the change in observed 
counts is 1 — 3%. Hence, it is safe to assume that the excess 
counts are solely caused by contamination by cluster mem¬ 
bers. Note, however, that the source redshift distribution is 
somewhat changed, as we do see intrinsically fainter galax¬ 
ies. We verified that the resulting change in mean redshift 
can be safely neglected. 

3.1 Contamination by cluster members 

To reduce contamination by cluster members, H12 used their 
limited color information to ide ntify and r e move galaxies on 
the red-sequence. As shown in iHoekstral (i2007f) . this does 
lower the contamination, but only by ~ 30% as many faint 
cluster members are blue. Furthermore H12 assumed that 
the excess number density of galaxies can be described as 
/contain oc r _1 , with the amplitude determined for each clus¬ 
ter. The analysis presented here differs from H12 in sev- 



radius [h^'Mpc] 


Figure 12. Plot of the ensemble averaged residual contamination 
as a function of distance to the cluster center. The top panel shows 
the results for sources with 22 < m r < 25 based on the Megacam 
data. The hatched regions indicate the 68% confidence intervals 
for the residuals for clusters with z < 0.25 (blue) and z > 0.25 
(red). The bottom panel shows the same, but now for CFH12k 
data and sources with 22 < m,R c < 25. For both data sets the 
residual contamination on the scales of interest (> O.S/iZTq 1 Mpc) 
is at most 2%. 


eral ways. Rather than applying a color cut, we restrict the 
magnitude range of the sources. Furthermore we use a more 
flexible model to quantify the radial dependence of the ex¬ 
cess counts. We also correct the ex cess source counts for 
t he o bscuration by cluster members dSimet fe Mandelbauml 
120141) . Finally, as described below, rather than considering 
the excess counts, we account for the weight provided by the 
uncertainty in the shape measurement. 

H12 used sources as bright as m r ~ 20, for which the 
level of contamination is high. This is demonstrated by Fig¬ 
ure EEU which shows the corresponding correction factor as 
a function of distance from the cluster for different bins of 
apparent magnitude. For the Megacam data the counts are 
normalised to the average number density at radii larger 
than r ma x = f/i^Mpc, i.e. we assume that the level of 
contamination can be ignored at those large radii. This is 
supported by comparison of the observed counts to those 
in blank fields. In addition, we used predic t ions b ased on 
the halo model described in ICacciato et al.l d2013l) to esti¬ 
mate the expected level of contamination due to neighbour¬ 
ing structures. In line with our comparison of the blank field 
galaxy counts, we find that the contribution from local struc¬ 
tures can indeed be ignored. The field-of-view of the CFH12k 
data is smaller and we estimate the background level using 
the number density at radii larger than r max = S/i^Mpc. 

The contamination is much higher for bright galaxies: 
such galaxies are rare in the field, but much more com¬ 
mon in clusters. For the brightest bin (20 < m r < 21) 
the cluster members outnumber source galaxies 3-to-l in 
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Figure 11. The ensemble averaged correction factor for contamination by cluster members as a function of distance for different bins of 
apparent magnitude. Panel (a) shows the results for clusters observed with Megacam, whereas panel (b) is for the CFH12k observations. 


the inner ~ l/i^Mpc. Such a large level of contamina¬ 
tion is difficult to model reliably and for this reason we 
decided to increase the bright limit of the source sample 
to m r = 22 from the typical value of m r = 20 used in H12. 
This leads to a reduction in excess counts that is comparable 
to excluding the galaxies on the red-sequence. Conveniently, 
the empirical correction for the multiplicative bias in the 
shape measurement also performs better for galaxies with 
22 < m r < 25. Furthermore, the lensing signal is higher 
for the fainter galaxies. The shapes of brighter galaxies are 
measured better and consequently given more weight in the 
lensing analysis (see Eqn.[5|. Rather than correcting for the 
excess counts, which effectively assumes that the weight is 
uniform, we compute the excess weight as a function of ra¬ 
dius. This is a minor correction, which increases the masses 
of the parametric NFW fits by ~ 2 — 3% (see > 14.11) . 

The other important change we make is that we allow 
the radial profile of the excess weight to vary from cluster to 
cluster by introducing a core. The simple 1/r profile used by 
H12 is not a good description for all clusters or magnitude 
bins. Investigation of the ensemble averaged residuals sug¬ 
gest that it leads to an overestimation of the contamination 
in the inner ~ 500/iyg 1 kpc and an underestimation by 4 — 5% 
at larger radii because the model attempts to compensate 
for the poor fit in the cluster cores. To describe the excess 
counts we now fit 


fc ontam (r) — 1 T U.Q ^ 


r + r c 


+ r c 


( 10 ) 


to each cluster, where we take r max = 4 /i^q 1 Mpc for the 
Megacam data and r max = Sh^Mpc for the CFH12k data. 
The core radius r c is a free parameter that we fit for each 
cluster separately. Figure [12] shows the ensemble averaged 
residual contamination for sources with 22 < m r /Rc < 2 5 
as a function of radius, suggesting that the systematic uncer¬ 
tainty for the ensemble of clusters is at most a few percent. 


Note that the residuals may be larger for individual clusters, 
resulting in increased scatter. However, the results presented 
in Figure [12] suggest that residual contamination will have a 
minimal impact on the normalization of scaling relations de¬ 
rived from CCCP measurements and we adopt a systematic 
uncertainty in the mass of 2% as a result of the imperfect 
correction for contamination by cluster members. 

The observed counts are biased low in the inner re¬ 
gions because the presence of bright cluster members af- 
fects our ability to detect a nd analyse sources. Although 
ISirnet fc Mandelbauml (I2014T ) showed that this is an impor¬ 
tant source of bias for the measurement of magnification, it 
may also lead to a small bias in our estimate of the dillution 
of the lensing signal. We simulated the impact of this and 
as described in more detail in Appendix El we find that the 
impact is indeed small, boosting the masses from the para¬ 
metric NFW fits by 1 — 2% (see > 14.ID . The aperture masses, 
which are discussed in > 14.21 are not affected, because they 
rely on estimates of the lensing signal at large radii where 
the density of cluster members is low. 


3.2 Photometric redshift catalog 

The weak lensing analysis of the initia l sample of 20 cl usters 
observed with the CFH12k camera in lHoekstral (120071 1 used 
the photometric redshift distributions derived for the HDF 
North and South using the availa ble deep multi-wavelength 
data (iFernandez-Soto et al.l[19991 1. However, the area cov¬ 
ered is small, leading to concerns whether the redshift dis¬ 
tributions are representative. For this reason H12 used th e 
photometric redshift distribution from lllbert et al.l (1200614 , 
which is based on the four CFHT Legacy Surve y Deep fields 
(each field covers one square degree). However. lllbert et al.l 
( 200614 derived photometric redshifts using observations in 
5 optical filters ( ugriz ). Although these data are very deep, 
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Figure 13. Left panel: The redshift distribution of galaxies with 20 < m r < 25 based on the COSMOS + UltraVISTA photometric 
redshift analysis. T he solid b lack h istogram indicates the results from our analysis, whereas the red dashed histogram corresponds to 
the results from lllbert et ahl ll2013l) . The inset shows the comparison of the photometric redshifts for the galaxies in common. The blue 
histogram shows the redshift distribution used in H12. Right panel: value of (8) as a function of m r for a cluster at 2 = 0.2. The points 
with error bars are our measurements, which agree well with th e red d ashed line ijllbcrt et alJbofii l for m r < 23. The four blue lines are 
the results for the CFHTLS Deep fields studied bv lllbert et alj J2006h and used in H12. 


the lack of near-infrared (NIR) data is a concern for high red¬ 
shift galaxies. Good quality NIR data are essential for this 
purpose because at 2 > 1.5 the Balmer and 4000A break fea¬ 
tures in galaxy spectral energy distributions, which are the 
strongest features for determining a photometric redshift, 
are redshifted into the NIR. 

For our analysi s we use data from th e Cosmic Evolution 
Survey ('COSMOS: IScoville et al.ll200'it )') which observed a 
single field covering 2 square degree with HST, and for which 
exten sive multi-w a veleng th and spectroscopic data are avail¬ 
able. Illbert et al.l (120091 4 present photometric redshifts for 
these data based on measurements in 30 bands (hereafter 
referred to as COSMQS- 3 0). T his redshift distribution was 
used in lAnnlegate et al.l l|2014l f to determine their “color- 
cut” masses. At the time of the analysis by Illbert et alJ 
d2009l4 . deep NIR data were not available. This situation has 
changed thanks to the UltraVISTA survey, an ESO public 
survey performin g deep imaging of the COSMOS field in 
5 NIR filters (see iMcCracken et al. 1 12011 for details). The 
UltraVISTA data are a significant upgrade to the available 
NIR imaging in the COSMOS field, and therefore allow for a 
marked improvement in the quality of photometric redshifts 
for galaxies at 2 > 1.5. Two public NIR-selected catalogs 
have been produced using the UltraVI STA data; one wher e 
galaxies were selected in the TG-band dMuzzin et al.|[2013lf . 
and one where objects w ere selected using a co-added \ 2 
image of the NIR bands (Illbert et al.l 120131 1. Both of these 
catalogs provide photometric redshifts and stellar masses for 
the galaxies. 

Galaxies have a wide range of optical-NIR colors, and 
therefore NIR-selected samples of galaxies are typically quite 
different from r-selected samples, particularly for the high- 


redshift end of the distribution. Consequently we cannot 
simply use the available photometric redshift catalogs, be¬ 
cause our source galaxies are selected from the deep CCCP 
r-band imaging. Furthermore, many of the sources that are 
of interest for the lensing analysis may be missing because 
the UltraVISTA data (K s = 23.9 AB) are 1.5 magnitudes 
shallower than the CCCP optical imaging. 

In order to construct a representative photometric red¬ 
shift distribution for the sources, we created a new r + - 
selected catalog of the COSMOS/UltraVISTA field using the 
Subaru imaging of the field from ICapak et al.l (120071 1. The 
Subaru r + imaging has good image quality (FWHM ~ OVS) 
and reaches a 5 a depth of ~ 26.5 AB. This is approximately 
a full magnitude deeper than the CCCP imaging, and there¬ 
fore provides a complete sample of galaxies to m r ~ 25.5 
with well-measured spectral energy distributions, a prereq¬ 
uisite for calculating photometric redshifts. The r + -selected 
catalog was constructed in th e identical ma n ner as the K s - 
selected catalog described in iMuzzin et al.l (l2013l l and we 
refer the reader to that paper for complete details of the cat¬ 
alog construction. In brief, the catalog consists of photom¬ 
etry in 29 photometric bands ranging f rom 0.15^tm — 24/im 
and inco rporates the availa ble GALEX (Martin et al. 2005), 
Subaru (ICapak et al.ll2007ll. UltraVISTA jMcCracken et al.l 
l2012lf . and Spitzer data ( Sanders et al.l 12007 1. Images were 
PSF-matched and photometry was performed in fixed 2'/l 
diameter apertures. 

Photometric redshifts for all galaxies were calculated 
using the EAZY photometric redshift code dBrammer et alj 
l2008h . which determines photometric redshifts using lin¬ 
ear combinations of multiple templates as well as a tem¬ 
plate error function to account for data/template mis- 
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Table 1. Basic information for the CCCP clusters 


(1) 

(2) 

name 

(3) 

z 

(4) 

RA 

(J2000.0) 

(5) 

DEC 

(J2000.0) 

(6) 

mag 

[h~g kpc] 

(7) 

</3> 

(8) 

(/3}ll3 

(9) 

(/^)used 

(10) 

5/3 

(11) 

(/3 2 )used 

1 

Abell 68 

0.255 

00 h 37 m 06.9 s 

+09 o 09'24" 

22-25 

0.553 

0.579 

0.566 

0.021 

0.375 

2 

Abell 209 

0.206 

01 h 31 m 52.5 s 

—13°36'40" 

22-25 

0.625 

0.649 

0.637 

0.019 

0.453 

3 

Abell 267 

0.230 

01 h 52 m 42.0 s 

+01°00'26" 

22-25 

0.586 

0.610 

0.598 

0.019 

0.409 

4 

Abell 370 

0.375 

02 h 39 m 52.7 s 

—01°34'18" 

22-25 

0.414 

0.442 

0.428 

0.026 

0.244 

5 

Abell 383 

0.187 

02 h 48 m 03.4 s 

—03°31'44" 

22-24.5 

0.636 

0.654 

0.645 

0.014 

0.462 

6 

Abell 963 

0.206 

10 h 17 m 03.8 s 

+39°02'51" 

22-25 

0.621 

0.644 

0.632 

0.018 

0.448 

7 

Abell 1689 

0.183 

13 h ll m 30.0 s 

—01°20'30" 

22-24.5 

0.647 

0.666 

0.656 

0.015 

0.475 

8 

Abell 1763 

0.223 

13 h 35 m 20.1 s 

+41° 00'04" 

22-25 

0.590 

0.613 

0.601 

0.018 

0.412 

9 

Abell 2218 

0.176 

16 h 35 m 48.8 a 

+66°12'51" 

22-24.5 

0.646 

0.662 

0.654 

0.012 

0.471 

10 

Abell 2219 

0.226 

16 h 40 m 19.9 s 

+46°42'41" 

22-25 

0.596 

0.621 

0.609 

0.020 

0.421 

11 

Abell 2390 

0.228 

21 h 53 m 36.8 s 

+17°41'44" 

22-25 

0.597 

0.624 

0.611 

0.021 

0.423 

12 

MS 0015.9+1609 

0.547 

00 h 18 m 33.5 s 

+16°26'16" 

22-25 

0.277 

0.304 

0.291 

0.025 

0.138 

13 

MS 0906.5+1110 

0.170 

09 h 09 m 12.6 s 

+10°58'28" 

22-25 

0.678 

0.700 

0.689 

0.016 

0.515 

14 

MS 1224.7+2007 

0.326 

12 h 27 m 13.5 s 

+19°50'56" 

22-25 

0.465 

0.492 

0.479 

0.024 

0.289 

15 

MS 1231.3+1542 

0.235 

12 h 33 m 55.4 s 

+15°25'58" 

22-25 

0.587 

0.614 

0.600 

0.021 

0.412 

16 

MS 1358.4+6245 

0.329 

13 h 59 m 50.6 s 

+62°31'05" 

22-25 

0.466 

0.494 

0.480 

0.025 

0.290 

17 

MS 1455.0+2232 

0.257 

14 h 57 m 15.1 s 

+22°20'35" 

22-25 

0.564 

0.594 

0.579 

0.024 

0.388 

18 

MS 1512.4+3647 

0.373 

15 h 14 m 22.5 s 

+36°36'21" 

22-25 

0.427 

0.458 

0.442 

0.027 

0.256 

19 

MS 1621.5+2640 

0.428 

16 h 23 m 35.5 s 

+26°34'14" 

22-25 

0.373 

0.404 

0.389 

0.027 

0.211 

20 

CL0024.0+1652 

0.390 

00 h 26 m 35.6 s 

+17°09'44" 

22-25 

0.393 

0.420 

0.407 

0.025 

0.226 

21 

Abell 115N 

0.197 

00 h 55 m 50.6 s 

+26°24'38" 

22-25 

0.645 

0.670 

0.658 

0.019 

0.478 


Abell 115S 

0.197 

00 h 56 m 00.3 s 

+26°20'33" 

22-25 

0.645 

0.670 

0.658 

0.019 

0.478 

22 

Abell 222 

0.213 

01 h 37 m 34.0 s 

—12 0 59'29" 

22-25 

0.620 

0.645 

0.633 

0.020 

0.449 

23 

Abell 223N 

0.207 

01 h 38 m 02.3 s 

-12°45'20" 

22-25 

0.629 

0.653 

0.641 

0.019 

0.459 


Abell 223S 

0.207 

01 h 37 m 56.0 s 

-12°49'10" 

22-25 

0.629 

0.653 

0.641 

0.019 

0.459 

24 

Abell 520 

0.199 

04 h 54 m 10.1 s 

+02°55'18" 

22-25 

0.642 

0.667 

0.655 

0.019 

0.475 

25 

Abell 521 

0.253 

04 h 54 m 06.9 s 

-10°13'25" 

22-25 

0.559 

0.583 

0.571 

0.021 

0.381 

26 

Abell 586 

0.171 

07 h 32 m 20.3 s 

+31°38'01" 

22-25 

0.668 

0.687 

0.678 

0.014 

0.501 

27 

Abell 611 

0.288 

08 h 00 m 56 8 s 

+36°03'24" 

22-25 

0.512 

0.536 

0.524 

0.022 

0.332 

28 

Abell 697 

0.282 

08 h 42 m 57.6 s 

+36°21'59" 

22-25 

0.532 

0.559 

0.545 

0.023 

0.354 

29 

Abell 851 

0.407 

09 h 42 m 57.5 s 

+46°58'50" 

22-25 

0.391 

0.418 

0.405 

0.026 

0.224 

30 

Abell 959 

0.286 

10 h 17 m 36.0 s 

+59°34'02" 

22-25 

0.528 

0.556 

0.542 

0.024 

0.350 

31 

Abell 1234 

0.166 

ll h 22 m 30.0 s 

+21°24'22" 

22-25 

0.693 

0.714 

0.703 

0.015 

0.534 

32 

Abell 1246 

0.190 

ll h 23 m 58.8 s 

+21°28'50" 

22-25 

0.651 

0.673 

0.662 

0.017 

0.483 

33 

Abell 1758 

0.279 

13 h 32 m 43.5 s 

+50°32'38" 

22-25 

0.539 

0.567 

0.553 

0.024 

0.361 

34 

Abell 1835 

0.253 

14 h 01 m 02.1 s 

+02°52'43" 

22-25 

0.562 

0.588 

0.575 

0.021 

0.385 

35 

Abell 1914 

0.171 

14 h 26 m 02.8 s 

+37° 49' 28" 

22-25 

0.685 

0.708 

0.697 

0.017 

0.525 

36 

Abell 1942 

0.224 

14 h 38 m 21.9 s 

+03°40'13" 

22-25 

0.607 

0.633 

0.620 

0.021 

0.434 

37 

Abell 2104 

0.153 

15 h 40 m 07.9 s 

—03°18'16" 

22-25 

0.707 

0.727 

0.717 

0.014 

0.552 

38 

Abell 2111 

0.229 

15 h 39 m 40.5 s 

+34°25'27" 

22-25 

0.599 

0.625 

0.612 

0.021 

0.425 

39 

Abell 2163 

0.203 

16 h 15 m 49.0 s 

—06°08'41" 

22-25 

0.619 

0.639 

0.629 

0.016 

0.445 

40 

Abell 2204 

0.152 

16 h 32 m 47.0 s 

+05°34'33" 

22-25 

0.708 

0.728 

0.718 

0.014 

0.554 

41 

Abell 2259 

0.164 

17 h 20 m 09.7 s 

+27° 40'08" 

22-25 

0.690 

0.711 

0.700 

0.015 

0.530 

42 

Abell 2261 

0.224 

17 h 22 m 27.2 s 

+32°07'58" 

22-25 

0.606 

0.632 

0.619 

0.020 

0.433 

43 

Abell 2537 

0.295 

23 h 08 m 22.2 s 

-02°11'32" 

22-25 

0.511 

0.537 

0.524 

0.023 

0.332 

44 

MS0440.5+0204 

0.190 

04 h 43 m 09.9 s 

+02°10'19" 

22-25 

0.646 

0.667 

0.656 

0.017 

0.476 

45 

MS0451.6-0305 

0.550 

04 h 54 m 10.8 s 

—03°00'51" 

22-25 

0.283 

0.307 

0.295 

0.024 

0.140 

46 

MS1008.1-1224 

0.301 

10 h 10 m 32.3 s 

—12°39'53" 

22-25 

0.504 

0.531 

0.517 

0.024 

0.326 

47 

RXJ1347.5-1145 

0.451 

13 h 47 m 30.1 s 

—11°45'09" 

22-25 

0.346 

0.371 

0.358 

0.024 

0.187 

48 

RXJ1524.6+0957 

0.516 

15 h 24 m 38.3 s 

+09° 57'43" 

22-25 

0.297 

0.321 

0.309 

0.024 

0.150 

49 

MACS J0717.5+3745 

0.548 

07 h i7 m 3Q 4 s 

+37° 45' 38" 

22-25 

0.269 

0.291 

0.280 

0.022 

0.131 

50 

MACS J0913.7+4056 

0.442 

09 h i3 m 45 5 s 

+40°56'29" 

22-25 

0.366 

0.393 

0.380 

0.026 

0.203 

51 

CIZA J1938+54 

0.260 

19 h 38 m 18.1 s 

+54°09'40" 

22-25 

0.550 

0.574 

0.562 

0.021 

0.371 

52 

3C295 

0.460 

14 h ll m 20.6 s 

+52°12'10" 

22-25 

0.343 

0.368 

0.356 

0.025 

0.185 


Column 2: cluster name; Column 3: cluster redshift; Column 4,5: right ascension and declination (J2000.0) of the adopted cluster center. 
In all but four cases (Abell 520, Abell 851, Abell 1758 and Abell 1914) we take this to be the position of the brightest cluster galaxy 
(BCG). Column 6: magnitude range used for the source galaxies. For clusters 1 — 20 this is the Rc filter and r' for the remaining clusters; 
Column 7: t he average value o f (3 = Di s /D s based on the photo-z analysis presented here; Column 8: the values for (3 measured as 
described in lllbert et alJ (|2013h : Column 9: the value for /3 we use to estimate masses, which is the average of the two measurements; 
Column 10: estimate for the systematic uncertainty in [3 as described in the text; Column 11: average value for (/3 2 ). 
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match. EAZY is well-tested and performs well amongst the 
best of publicly-availab le photometric redshift codes (e.g. 
IHildebrandt et al.l 120101 1. The photometric redshifts were 
further refined by determining small offsets to the photo¬ 
metric zeropoints using the ~ 5000 spectroscopic redshifts 
available in the COSMOS field from the zCOSMOS-lOk 
sample dLillv et al.ll2009l ). The process is iterative, and the 
final photometric catalog contains photometric redshifts ac¬ 
curate to Az/(1 + z) = 0.01, with a catastrophic outlier 
fraction of ~ 1%. This estimate of the accuracy is based on 
the zCOSMOS spectroscopic redshifts, which are primarily 
bright galaxies at z < 1.5. At z > 1.5 it is more difficult to 
assess how accurate the photometric redshifts are, due to the 
lack of a spectroscopic calibration sample. Small numbers of 
spectroscopic redshifts are available from various NIR spec¬ 
troscopic surveys, and those suggest that the accuracy at 
z > 1.5 for bright galaxies is only slightly worse, of order 
a few percent. We have made the full r + -selected catalog 
publicly available on the As-selected catalog websitcQ Also 
included is a simplified version of that catalog that can be 
used for a quick calculation of photometric redshift distri¬ 
butions for future lensing analyses. 

The left panel in Figure [13] shows the resulting red¬ 
shift distribution for galaxies with 20 < m r < 25 us¬ 
ing our r-band selected photometric redshift catalog. The 
red dashed histogram shows the redshift distribution of 
the galaxie s matc hed to the NIR-selected catalog from 
lllbert et al.l (l2013h . which makes use of the ~ 3000 un¬ 
published zCOSMOS-deep spectroscopic redshifts. For ref¬ 
erence, the blue histogram shows the correspon ding redshift 
distribution for the CFHTLS Deep fields from lllbert et alj 
(120061 1 . which was used by H12. The inset panel shows a 
direct comparis on of the ph o tomet ric redshifts derived here 
and those from lllbert et all (120131 ) who used a different al¬ 
gorithm and measured the photometry independently. The 
overall agreement is remarkably good for the galaxies in 
common; only for m r > 24 are some of the galax i es ass igned 
a high redshift in the catalog from lllbert et alj d2013lf and 
a low redshift by EAZY. These represent 2% of the total 
sample of sources and about 20% of the galaxies for which 
lllbert et al.l d2013h find z > 2. The excellent agreement for 
most of the galaxies is a demonstration of the quality of the 
data, both in terms of depth and wavelength coverage. 

The impact of the differences in source redshift distribu¬ 
tions on the cluster mass estimates is quantified in the right 
panel of Figure fHfl which shows the lensing efficiency (/3) as 
a function of apparent magnitude. Our results are indicated 
by the black points, whereas the red dashed curve corre - 
sponds to the redshift distribution from lllbert et all d2013h . 
For m r < 23 the agreement is ver y good, whe r eas th e higher 
number of z > 2.5 galaxies in the lllbert et al.l (l2013l i catalog 
results in a higher value for (/3), and thus a lower mass. We 
list (/9) for the sources in Table [T| Column 7 lists the values 
for our analysis of the COSMOS and UltraVIS TA data, and 
Colum n 8 lists the results using the results from lllbert et all 
d2013lf . To determine cluster masses we use {/3) US ed provided 
in Column 9, which is the average of the values obtained 
for the two redshift distributions. The value for (/3 2 ) US ed, 
which is a measure of the width of the distribution (see 


5 http://www.strw.leidenuniv.nl/galEuxyevolution/ULTRAVISTA/ 


iHoekstra et al.l l200d : iHoekstral 120071 . for details) is also an 
average of the two estimate^]. These numbers include a size 
cut similar to the one applied to our lensing data. The red 
line in Figure [14] shows how (/?) depends on the observed 
half-light radius; we End that the dependence on size is very 
small, and we therefore conclude that the size cuts do not 
introduce a significant bias. 

The significant differences in redshift distribution 
demonstrate that the lack of reliable photometric redshift 
estimates remains a key source of error. The unique range 
in wavelength and quality of the COSMOS+UltraVISTA 
data are a major step forward, but without complete spec¬ 
troscopic coverage, the uncertainty at the highest redshifts 
remains. Furthermore, cosmi c variance can still be impor¬ 
tant for a single field (see e.g. IHoekstra et aljfeoill . for esti¬ 
mates). Assigning a systematic uncertainty remains difficult, 
but we use the difference be t ween our photometric redshifts 
and those from lllbert et all (120131 / (indicated as 113 in Ta¬ 
ble [TJ as an estimate for the systematic uncertainty. 

The blue lines in the right pan el of Figure [T51 show {/3) 
for the four CFHTLS Deep fields (lllbert et al.l 20061 1. The 
main difference occurs at m r ~ 21.5, although the average 
is also lower than the new estimates for faint galaxies. We 
use the variation between the fields to estimate the contribu¬ 
tion of cosmic variance for our redshift distribution. This is 
conservative, because the COSMOS field itself is larger than 
each of the CFHTLS Deep fields. The systematic uncertainty 
S/3 listed in Table [l] is the sum of the dispersion measured 
from the four fields and |/3ii3 — /3|/2. These amount to ~ 3% 
for clusters at z = 0.2, increasing to ~ 8% for z = 0.55. 

Variation in the actual source redshift distribution be¬ 
hind a cluster leads to an increase in the st atistical uncer¬ 
tainty . The impact of this was studied by IHoekstra et al.l 
(1201 ll) using simulations. Their findings suggest that the 
lack of redshift information for the individual sources in¬ 
creases the statistical uncertainty in the mass by ~ 3% for 
a cluster at z = 0.2 and by ~ 10% at z = 0.6. 


4 UPDATED MASS ESTIMATES 

We use our new insights in the shape measurement bias and 
the source redshift distribution to update the mass estimates 
provided by H12. The only change we made to the original 
shape measurements is a corrected estimate for a e , the un¬ 
certainty in the polarizatiorQ because we use this quantity in 
our estimate of the multiplicative bias. This correction also 
affects the weight defined by Eqn. [5] although the impact 
is very minor. The star-galaxy separation is somewhat more 
conservative: previously we included faint galaxies with sizes 
comparable to the PSF in the lensing analysis, although they 
were severely downweighted in practice. As the correction 
scheme requires sizes larger than the PSF, and the fact that 
these objects do not contribute much to the average signal, 


6 The values listed here are corrected for an error in the calcu¬ 
lation of (/3 2 ) that reduced the senstitivity to the convergence in 
our previous work. 

7 The older version of the code included an incorrect treatment 
of the Poisson noise. 
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Figure 14. The ensemble averaged lensing signal (rg/)3) as a 
function of observed source half-light radius . The value for the 
Einstein radius r e was obtained by fitting a SIS model to the 
lensing signal at radii 0.75 — 4//7 q' Mpc. The solid points, which 
should not depend on source size, indicate the results when we 
apply our correction for the multiplicative bias, whereas the open 
points are the uncorrected values. The lensing signal is biased 
low for galaxies with ry, > 5, similar to what we observed in 
our simulated data. The red line indicates the estimate of (ft) as 
a function of source size (the values are indicated on the right¬ 
most vertical axis), which does not vary significantly with object 
size, suggesting that the adopted redshift distribution is not very 
sensitive to the size cuts applied. 


we now only select objects larger than the PSF in the ob¬ 
ject catalog. The star selection itself was not changed, and 
consequently we use the same PSF model parameters. 

We apply our empirical correction given by Eqn. [9] to 
the object catalogs using the observed values of and v and 
recompute the tangential shear profile as a function of ra¬ 
dius, taking the cluster centers used in H12 (Table 1 of that 
paper). Based on our image simulations we also apply a size 
cut rh < 5 pixels. The solid points in Figure [14] show the 
lensing signal, averaged over the full ensemble of clusters, as 
a function of the observed half-light radius of the sources. 
The signal is quantified by the Einstein radius te obtained 
from a fit to the measurements between 0.5 — 4/iyQ 1 Mpc. We 
correct each bin for the mean ft, which does not vary signifi¬ 
cantly with source size (indicated by the red line), suggesting 
that our results are not sensitive to the size cuts. The re¬ 
sulting value of r_g/ ft should not depend on the source size. 

The observations support our finding from the image 
simulations that the signal is biased low for galaxies with 
m > 5 pixels. For reference we also show the results if we 
do not apply our empirical correction for multiplicative bias 
(open points). The corrected results are consistent with a 
constant signal, but the uncorrected measurements are bias 
low for the smallest sources. 

We update the correction for G alactic extinction us¬ 
ing the ISchlaflv fc Finkbeined (1201 ll ) recalibration of the 


ISchlegel et al.l Jl998l ) infrared-based dust map. This reduces 
the correction for clusters in highly extinguished regions. To 
ensure a more robust correction for the contamination by 
cluster members we do not apply a color cut, but instead 
limit the source sample to galaxies with 22 < m r /R c < 25. 
In the case of the CFF112k data, H12 used galaxies fainter 
than Rc = 25, for which the redshift distribution is not well 
constrained. We now limit the analysis to galaxies brighter 
than Rc = 25. 


4.1 Parametric mass models 


One approach to infer masses is to fit parametrized mod¬ 
els to the lensing signal. The most commonly used profile is 
th e Navarro-Frenk-Whi te (NFW) fitting function proposed 
by I Navarro et all (1 19971 ). which is a good description of the 
average density profiles of halos in numerical simulations 
of structure formation in cold dark matter dominated uni¬ 
verses. It also describes th e stacked lensing si g nal for ensem¬ 
bles of clusters well (e.g. lOkabe et al.l 120131 : lUmetsu et alJ 
120141) . The NFW profile is characterized by two parameters. 
We use the mass of the halo and the concentration c, al¬ 
though we do not fit for these parameters simultaneously: 
we use the fact that simulations show that the mass and con- 
centraton are correlated. However, as the concentration de¬ 
pends on the formation redshift of the halo, this relation de¬ 
pends on cosmology. H12 used the results from lDuffv et all 
|2008j) , which are based on the cosmological parameters from 
the five-year Wilki nson Microwave Anis otropy Probe obser¬ 
vations (WMAP5: iKomatsu et aI1l2009l ). These have since 
been superseded by the measurements o f the Planck satel- 
lite ( Planck Collaboration et ~ahl l2014bl ). iDutton fe Maccibl 
toil present fitting functions for the mass-concentration 
relation for this cosmology, which we use when we estimate 
clus ter masses (cf. Table l2ll. 

iBecker fe Kravtsovl (|201ll ) have shown that fitting an 
NFW model to the observed lensing signal can lead to mass 
estimates that are biased low when measurements at large 
radii are included. For this reason we restrict the fit to 
0.5 — 2 hjQ Mpc from the cluster, where biases should be 
negligible. The resulting massed, for different overdensities 
A are presented in columns 9 — 11 in Table [2] The statisti¬ 
cal uncer tainties on the m easurements are estimated as de¬ 
scribed in lHoekstral (120071 ') and H12. The uncertainties in the 
mass estimates include the contribution from distant large - 
scale structure (lHoekstrall200ll . 120031 : iHoekstra et a.l .11201 ll ) . 
For reference with other studies we also present the velocity 
dispersion corresponding to the best-fit singular isothermal 
sphere (SIS). 

For reference we note that i f we had u sed the mass- 
concentration relation from lDuffv et al.l d2008l) . which yields 
concentrations that are ~ 20% lower compared to the val¬ 
ues used here, our masses would change as follows: M2500 
decreases on average by 7%, while M500 and M200 increase 
by 5% and 9%, respectively. The relative change in mass 


8 Ma is the mass enclosed within a radius where the mean den¬ 
sity of the halo is A times the critical density at the redshift of 
the cluster; the virial mass is defined relative to the background 
density. See lHoekstral d2007l ) for more details about our choice of 
definition. 
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Table 2. Weak lensing mass estimates for the CCCP sample 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ii) 

name 

<7 

MET 

met 

7*2500 

M 2 500 

7*500 

M 500 

m nfw 

vir 

m nfw 

1V1 2500 

mZ w 

Abell 68 

H17±n 

4.5 ±0.5 

8.3 ±1.3 

552 

0 1+0.4 
°' L -0A 

1380 

9.7++° 

42.91!;! 

0 0+0.6 

0 -0.6 

7-sll:! 

Abell 209 

970+™ 

3.2 ±0.5 

7.3 ±1.2 

490 

2 l+°- 4 
Z ’ 1 -0.4 

1230 

6-5±i1 

«-4l!:2 

9 0 + 0.5 
Z ' U -0.6 

4 9”*" 1 ' 3 

Abell 267 

1006+ gj 

3.9 ±0.5 

7.2 ±1.4 

525 

2 6” 1 ” 0 ' 4 
z '°-0.4 

1203 

6-3±i;l 

7 9"*" 2 ' 4 

1 *- 2.3 

1-918:8 

4-7!l;l 

Abell 370 

1489179 

6.9 ±0.7 

15.3 ± 1.6 

638 

r r + 0.8 
°'°-0.8 

1637 

18.51!;! 

30.41!;! 

6-511:! 

17.5!!;! 

Abell 383 


2.5 ±0.6 

6.1 ± 1.3 

430 

1 4+0-6 
1 ' 4 -0.5 

1217 

6 + 2 ! 

5-si!;! 

1 4+0-6 
1 '^-0.6 

3 4+ 1 - 5 

Abell 963 

11061 ™ 

3.4 ±0.5 

6.6 ± 1.4 

506 

9 q+0.5 
^-0.4 

1185 

5-8+5 

1 9 q+3.0 
lz '°-3.0 

2.918+ 

7 1+ 1 - 7 
* " L —1.7 

Abell 1689 

1429lg 

7.0 ±0.5 

13.2 ± 1.4 

702 

5 9 + 0-7 
°' y -0.7 

1571 

13.31!+ 

30.91!! 

e-ell:! 

17 . 3 !!;! 

Abell 1763 

12291™ 

4.9 ± 0.6 

10.0 ± 1.4 

604 

3 q+ 0-6 
°' y - 0.6 

1511 

-I 9 q + 3.2 

i6.9l!;! 

q q+ 0.8 
^' 8 - 0.8 

9 7 + 2.0 

1 - 2.0 

Abell 2218 

11811™ 

5.0 ±0.6 

8.7 ±1.4 

630 

4 2 +0 ' 7 
^-0.7 

1379 

8 9 +2,2 
°' y - 2.1 

16.41!* 

q 7+0.9 

^•'-o.s 

9-41®;! 

Abell 2219 

10411™ 

4.1 ±0.5 

9.5 ± 1.4 

552 

3 Q+ 0- 5 
4- u -o.5 

1408 

10 . 0 l™ 

11 3 +2 ' 5 
11 , 0 -2.5 

9 7 +O .6 

Z ‘' — 0.6 

6 - 611 ;! 

Abell 2390 

133llg4 

4.9 ±0.5 

9.9 ± 1.3 

602 

3 g+ 0- 5 
°' y -0.5 

1351 

S-Sli'J 

9 q 1 +3.8 
^°'-3.6 

5 1 +O .8 
°' 1 - 0.8 

13.2l!;l 

MS 0015.9+1609 

145611™ 

6.8 ± 0.6 

17.9 ± 1.9 

601 

°'- 0.8 

1617 

2 i- 8 l!;! 

28.91™ 

6 2”*” 1 ' 5 
°' Z — 1.5 

16 9 +41 

lu ' y -4.0 

MS 0906.5+1110 

10771™ 

3.8 ±0.5 

8.7 ± 1.2 

549 

2 . 8«; 5 5 

1423 

9-7l T 6 

12 . 61 !° 

2 9 " 1 ” 0 ' 7 
- 0.6 

7 3 ”*” 1 ' 7 

1 *°-l.6 

MS 1224.7+2007 

860lj3g 

1.8 ± 0.6 

2.4 ± 1.8 

345 

0 - 8 io .7 

782 

1 9 + 1,0 
± ' b, -0.9 

4.9l!;? 

-1 9 + 0.6 

1 ‘ z — 0.5 

3 O "*" 1 ' 4 
°' u -1.3 

MS 1231.3+1542 

5901™® 

1.0 ±0.5 

0.4 ± 1.3 

344 

0 7 + 0 - 2 
u -' - 0.2 

565 

o-e± 8 :S 

i-9ll:! 

0 c+0.3 
U '°-0.3 

x 4 + 0.8 
1 ' x —0.7 

MS 1358.4+6245 

11671™ 

4.3 ±0.5 

8.7 ± 1.5 

529 

3 0+0-5 
°' u -0.5 

1291 

8 - 611:8 

i3-4l!:! 

0 i+0.7 
± —0.7 

7 0+1-9 

7 ' 8 — 1.8 

MS 1455.0+2232 

H31+ 63 

3.7 ±0.4 

7.2 ± 1.2 

510 

9 c+0.4 
z '°-0.4 

1158 

5 -7li;! 

13.21!+ 

q -| + 0.6 
°'- l - 0.6 

7 7 + 1-4 

‘ *' -1.5 

MS 1512.4+3647 

700 + 111 
'°°-130 

1.5 ±0.6 

4.5 ± 1.4 

282 

0 tr+0.3 
U ‘ D -0.3 

853 

Z '°- 1.6 

3.911:1 

1 o +0 - 4 

±>u -0.4 

2 4 4 " 1 - 1 
z '^-i.i 

MS 1621.5+2640 

1300lgg 

4.9 ± 0.6 

11.4 ± 1.7 

543 

0 0 + 0.8 
o,o - 0.8 

1286 

9-511:8 

19.111! 

4 o+O.9 

^'°-0.9 

11 2”*” 2 ' 5 
11 ' z -2.4 

CL 0024.0+1652 

13111& 

5.6 ±0.6 

11.4 ± 1.7 

571 

4 0 + 0-6 
^• u - 0.6 

1333 

10 . 21 !+ 

24 4 + 5,7 
z ^-5.4 

r 0 + 1.3 
°'°-1.2 

14 l - *" 3,3 
1 ^' ± -3.1 

Abell 115N 

oqo+89 

o>JO_99 

1.4 ±0.4 

5.3 ±1.1 

283 

0 4 +°- 3 
u ‘^-0.4 

1098 

4.6ll° 

5 9 + 20 

1 c+0.5 
1 , 0 -0.5 

q c + 1-2 
o.°_i 2 

Abell 115S 

859l®j 

2.6 ±0.4 

5.9 ± 1.2 

416 

-1 9 +O .4 
1 ,Z -0.5 

1127 

s-oll;! 

7 O -1 ” 2 ' 0 

1 ' u - 2.0 

1 7 + 0.5 
± * -0.5 

4.111;! 

Abell 222 

9161^ 

2.9 ±0.5 

6.9 ± 1.4 

450 

1 6~*~ 0 ‘ 6 
1 -°- 0.8 

1174 

5-711;! 

5 4 + 2.1 

0 '^-2.1 

l 0+0.5 
i '°-0.5 

3.811;! 

Abell 223N 

989l*° 

3.0 ±0.5 

7.5 ± 1.3 

463 

1 7+0.5 
1 -'-0.6 

1236 

e-ell;! 

« 9+ 2 - 5 

°' y -2.5 

9 -1 + 0.6 
z ' 1 - 0.6 

tr 0+1-5 
°' Z -1.5 

Abell 223S 

92319° 

3.0 ±0.5 

8.3 ±1.1 

466 

1 S+°- 5 
L8 —0.5 

1370 

y ' u -1.5 

7 o+ 2.6 

1 '°-2.4 

1 g+0-6 
x ' y -0.6 

4-6ll;| 

Abell 520 

1144l« 

3.6 ±0.4 

7.3 ± 1.1 

526 

9 c+0.5 
z '°-0.4 

1208 

e-ill:? 

1 15 o+3.0 
lo '°-3.0 

0 c+0-7 
°'°-0.7 

8.811;! 

Abell 521 

9441°^ 

3.1 ±0.5 

9.0 ± 1.4 

448 

1 7+0-7 
J -'- 0.8 

1335 

8 - 8 l ?;8 

10 7+ 30 
' —3.0 

2 c i +0 - 7 
Z '°-0.7 

6.311;! 

Abell 586 

8041^2 

2.5 ±0.6 

6.1 ± 1.6 

441 

1 4 + 0-5 
± ’^ t — 0.4 

1221 

6 l" 1 " 2,6 
u ' 1 —2.6 

4 6" 1 " 2 ' 1 
^' u -2.0 

1 9 +O .5 
4 ‘ Z — 0.5 

2.811! 

Abell 611 

995l?o3 

3.7 ±0.5 

9.0 ± 1.4 

502 

2 4 +0 - 5 

Z '^-0.5 

1236 

7 9 +I .5 
' * z — 1.4 

9 4+ 2 - 8 
y '^-2.8 

9 9 +O .7 
Z ' Z — 0.7 

5-511;! 

Abell 697 

11461™ 

4.6 ± 0.5 

10.5 ± 1.4 

565 

3 4+0-6 
0.6 

1431 

11 . 2 II® 

14 l+ 31 
3.1 

q q+0.7 

0.7 

8 . 211 ;! 

Abell 851 

13281°! 

5.4 ±0.6 

12.2 ± 1.6 

553 

q 7+0.5 

1362 

11 1+ 2 - 2 
2.1 

21 4+ 5 - 3 
Z 1 '^-4.9 

4 7+ 1 - 2 

-1.1 

12-51!;! 

Abell 959 

12571™ 

5.0 ±0.5 

10.8 ± 1.4 

596 

4 0+0-6 
^• U -0.6 

1343 

9 2 +1 ' 6 
y * z — 1.6 

19 7+3-8 
-3.6 

4 4+0-8 
^- 0.8 

11 4“*" 22 

Abell 1234 

9691™ 

2.5 ±0.5 

4.4 ± 1.3 

447 

1 c+0.3 
1 -°-0.3 

989 

q 9+1-2 
°' z -1.0 

7-el!:! 

1 9+ 0 - 6 

4 5 "*" 1 ' 3 
^. 0-12 

Abell 1246 

9211 ll 

2.7 ±0.4 

5.6 ±1.1 

440 

1 c+0.3 
± '°-0.4 

1089 

4-411;! 

s-7l!:! 

9 -1 + 0.6 

Z ' 1 -0.5 

51+ 1 - 4 
°' 1 —1.3 

Abell 1758 

12781“ 

5.5 ±0.5 

12.1 ± 1.4 

651 

c 9+0- 7 
°' Z -0.7 

1507 

12.911;! 

i9.4l!:! 

4 q+0- 7 

11 2”*” 1 ' 9 
11 * z — 1.9 

Abell 1835 

12951“ 

5.3 ±0.5 

10.7 ± 1.3 

618 

4 3 + 0-5 
^- 0.5 

1398 

io.oll;! 

19.91 ll 

4 4+0-8 
^'^- 0.8 

ii-4l!:! 

Abell 1914 

10981™ 

3.7 ±0.5 

7.8 ± 1.2 

531 

9 tr+0.4 
z '°-0.4 

1293 

7 q+1-3 
‘ '°-1.3 

13.51!+ 

q 1 +0.6 
°“ L — 0.6 

7-811:1 

Abell 1942 

10801™ 

3.8 ±0.6 

7.5 ± 1.3 

531 

9 7 +O .6 

Z '' —0.5 

1212 

6-411;! 

i3.5l!:! 

q 1 +0.6 
°“ L — 0.6 

7 0+1-6 

7 * 8 —1.5 

Abell 2104 

11351™ 

4.2 ±0.5 

10.2 ± 1.2 

596 

q c+0.6 
°' D -0.6 

1426 

9-611.1 

i5-7l!:! 

0 0+0.8 
' 3 '°-0.8 

9 0 + 2-0 
y ' u -1.9 

Abell 2111 

9961™ 

3.9 ±0.5 

6.6 ± 1.5 

528 

2 6” 1 ” 0 ' 5 
Z '°-0.5 

1170 

r 7+1-9 
- 1.8 

9-4l!:! 

9 q+0.6 

z ’°—0.5 

c c + 1-4 
°'°-1.3 

Abell 2163 

1188179 

4.4 ±0.4 

9.4 ± 1.2 

574 

q q+0.4 

—0.4 

1466 

n-ol !:8 

17 4+ 3-8 

0 g+0.9 
-0.8 

10 . 0 !!;! 

Abell 2204 

12291“ 

4.8 ±0.5 

10.8 ± 1.2 

631 

4 o+O - 5 
^' z - 0.6 

1491 

u-oll:S 

19.91!;! 

4 4+°- 7 
^'^-0.7 

11-311 f 

Abell 2259 

Q32 + 89 
yo ^-98 

2.4 ±0.6 

5.6 ± 1.2 

427 

-1 q + 0.5 
- l ' o -0.4 

1113 

4-ell;! 

7 9+ 2 -6 
' ' y -2.4 

1 9+0-6 
x ' y -0.6 

4-611:1 

Abell 2261 

13071“ 

6.0 ±0.5 

14.1 ± 1.4 

682 

c 7 +O .6 

O' ' —0.6 

1663 

16.411;! 

24-41!+ 

c: q+0. 9 

-0.9 

i3.9l!! 

Abell 2537 

12851™ 

5.4 ±0.6 

10.0 ± 1.4 

599 

4 1 + 0-6 
^• X - 0.6 

1312 

8 -7ll;! 

20 . 9 !!;! 

4 6" 1 ” 0 ' 9 
^'°-0.9 

12 l - * -2,3 
lz,± - 2.2 

MS 0440.5+0204 

78011™ 

2.9 ±0.6 

2.8 ± 1.3 

468 

1 q+0-6 

1 ' 8 -0.5 

896 

2 5 + 0 - 7 
z '° —0.7 

3-511;! 

0 g+ 0- 5 
u ' y -0.5 

2 -i!l:l 

MS 0451.6-0305 

1302l™| 

4.4 ±0.6 

8.4 ±2.1 

466 

2 6~*~°' 6 
Z '°-0.5 

1082 

D '°-1.9 

i7-4l!;! 

3 9 ”*” 1 ' 3 
°' y - 1.2 

10.3!!;! 

MS 1008.1-1224 

12181 ™ 

4.1 ±0.4 

8.2 ± 1.4 

520 

9 ' 7 + 0.4 

z ' '-0.4 

1176 

°'°- 1.2 

i6.3l!;! 

q 7+0.8 

1 - 0.8 

9 5 +2 '° 
y '°-1.9 

RX J1347.5-1145 

1358l“ 2 

5.2 ±0.7 

10.1 ± 2.2 

547 

q o + 0.9 

^ -0.9 

1323 

40-71!;! 

19.9!!;! 

4 4+ 1 - 2 

11 7 -1 " 3,0 
' —3.0 

RX J1524.6+0957 

8391™! 

2.1 ±0.7 

6.5 ±2.1 

249 

0-4 ±? 2 3 2 

977 

4 6 +2 ’ 2 
^' u - 2.0 

6 - 2 l!;? 

i-5!l:! 

q 0 + 2.5 

4 - 0—2 2 

MACS J0717.5+3745 

16171™! 

6.4 ±0.8 

19.3 ±2.3 

612 

s-8ti:5 

1489 

17 l + 3 - 2 

1 ' 3.1 

38.4!!;! 

8 .o!l;| 

99 q + 5.2 

^•°-5.0 

MACS J0913.7+4056 

919+ 143 
^ -168 

3.1 ±0.8 

5.3 ± 1.8 

397 

1 4+0 - 9 
-0.7 

945 

3 9 + 1,3 

1.2 

6-511:8 

1 6" 1 ” 0 ' 8 
i '°-0.7 

0 g + 2.0 
°' y — 1.8 

CIZA J1938+54 

11861 ™ 

5.3 ±0.6 

11.3 ± 1.5 

601 

4 0 + 0-6 
^' u - 0.6 

1573 

14 4 + 2 - 4 
± ^' -2.4 

17.211;! 

0 g+0.9 

°' y -0.9 

9 9 -1 " 2 ’ 4 

3C295 

107611™ 

4.6 ±0.7 

8.1 ± 1.9 

501 

2 9 +0,7 
z ' y - 0.6 

1101 

6 . 21 I! 

12 . 0 !!° 

9 o+O.9 

0.9 

7 1 +2.3 
' 1 — 2.2 


Column 1: cluster name; Column 2: line-of-sight velocity dispersion (in units of km/s) of the best fit SIS model; Columns 3 & 4: projected 
mass within an aperture of 0.5 /i^q 1 and lh^Mpc, resp.; Columns 5 & 7: r/\ (in units of h^kpc) determined using aperture masses; 
Columns 6 & 8: deprojected aperture masses within t*a; Columns 9-11: masses from best fit NFW model. All masses are listed in units 
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does not depend significantly on the cluster redshift. How¬ 
ever, for a direct comparison with the existing li ter ature , 
we present mass measurements using the lDuffv et alJ ( 2008 1 
mass-concentration relation for the WMAP5 cosmology in 
Appendix m 

4-1-1 Comparison to other weak tensing studies 

Several studies have determined weak lensing masses for 
large samples of clusters using observations with the Sub¬ 
aru telescope. The most relevant for the comparison with 
CCCP is the Weighing the Gian ts (WtG) project, described 
in Ivon der Linden et al.l (l2014al l. which targeted 51 massive 
clusters. For a subset of the clusters WtG determined p ho- 
tometric redshifts for the sources fsee lKellv et al.ll2014 . for 
details). However, these are predominantly the high redshift 
systems where the overlap with CCCP is limited. For this 
reason we c ompare with the “color -cut” masses, which are 
presented in I Applegate et all (l2014l l. 

A closer inspe c tion o f th e sample studi e d by 
Ivon der Linden et al.l (l2014al ) and lApplegate et al.l (120141 ) 
shows that they associated MS0906.5+1110 with the clus¬ 
ter Abell 750 which is lo cated on ly 3’ away in projection. 
However, as discussed bv lRines et alJ (I2013T) the latter is a 
different cluster, which is cle arly separated in redsh i ft. The 
location of A750 provided in Ivon der Linden et al.l d2014al l 
is in fact that of MS0906.5+1110, and we therefore include 
this cluster in the comparison. Abell 1758 is a merging clus¬ 
ter and therefore H12 considered the Eastern and West - 
ern component separately (also see iRagozzine et al.l 120121) . 
However, other studies consider this a single cluster and we 
th erefore decided to pro v ide res ults for the location listed 
by Ivon der Linden et al.l d2014al l , who refer to this cluster 
as Abell 1758N. As a result we have 18 clusters in common 
with WtG. 

To compare the results for these clusters, we follow 
lApplegate et al.l d2014l l and fit an NFW model with a con¬ 
centration C200 = 4 to the tangential distortion within 
0.75 — S/i^Mpc and compute the mass within a sphere of 
radius 1.S/i-/, 1 Mpc. The results are presented in Figure [To] 
We find that the WtG masses are somewhat larger: the 
dashed line indicates the best fit linear relation MwtG = 
(1.082 ± 0.038) Mcccp • Repeating the comparison using the 
results from H12 yields MwtG = (1.263+0. 048 )Mhi 2 - Hence 
the analysis presented here reduces the discrepancy con¬ 
siderably. We note that differences in the fitting procedure 
can lead to additional uncertaintly, and it is therefore not 
clear whether the differ ence is significant . Furth ermore, the 
“color-cut” masses from lApplegate et aid d2014l ) are derived 
using the photometric redshift catalog from lllbert et aid 
d2009h . which are based on the original COSMOS-30 data. 
Using this redshift distribution we find MwtG = (1.063 ± 
0 . 038 )Mcccp. Interestingly, when we compute deprojected 
aperture masses within a radius of l/i^Mpc (in this case 
adopting C200 = 4; see 1)4.21 for details), the agreemen t with 
the corresponding masses from lApplegate et alJ (120141 ) is ex¬ 
cellent: MwtG — (1.018 ± 0.036)Mqq C p . 

The NFW model is fit to relatively small radii, where 
the contamination by cluster members is large (although the 
inner 750 /i/q 1 kpc are not used): if we omit the correction for 
the contamination of cluster members our masses decrease, 
as expected, and MwtG ~ 1-28 x Mcccp- Although the 



M cccp (<1.5 h 70 ‘ Mpc) [10“ h 7 - 0 > M 0 ] 


Figure 15. Compa rison with the mass estimates from 
lApplegate et all (120141 ). The CCCP masses are computed from 
the best fit NFW model to the lensing measurements within 
0.75 — 3/iyQ 1 Mpc, adopting a concentration C200 = 4 . This matches 
the procedure described in lApplegate et al.l ll2014ll . although our 
source redshift distribution is somewhat different, as explained in 
the text. The dotted line indicates the line of equality, whereas 
the dashed line is the best fit, which has a slope 1.08. 


correction is substantial, Figure Q3 suggests that the bias 
after correction should be < 2%. 

We investigated this further by restricting the fit to 
small (r in = 0.75— l.Shjg Mpc) and large (r out = 1.5 —3 /i/q 1 
Mpc) radii. If the contamination correction is adequate, 
the resulting average masses should agree, whereas a ratio 
M ou t/Mi n > 1 would imply residual contamination by clus¬ 
ter members. For the 18 clusters in common with WtG we 
find M ou t/Mi n = 1.05 ± 0.05, suggesting that the correction 
has worked well (the ratio is 1.16 + 0.05 if we do not cor¬ 
rect for contamination). For the full CCCP sample we find 
M ou t/Min - 1.00 + 0.03 

lUmetsu et al.l (120141 ') present results for 20 clusters stud¬ 
ied as part of the Cluster Lensing and Supernova survey with 
Hubble (CLASH). Of these 17 clusters overlap with WtG, 
but only 6 overlap with CCCP. For the six clusters we have 
in common, the CLASH masses are 12 + 5% hig her than the 
CCCP results. However, we note that lUmetsu et al.l d2014l ) 
use a different fit range, while leaving the concentration a 
free parameter. Although they find a best fit concentration 
of ~ 4 when they stack the clusters in their sample, cluster- 
to-cluster variation complicates a more direct comparison. 
As discussed below, we analysed the CLASH data using our 
pipeline and find better agreement with our masses derived 
from CFHT observations. 

Another large study that does overlap considerably with 
CCCP is the Local Cluster Substructure Survey (LoCuSS). 
Results for 3 0 clusters, of which 1 3 overlap wit h CCC P are 
presented in lOkabe et al.l (l2010l) . lOkabe et al.l ll 20131) sug¬ 
gest that a revised analysis leads to higher masses, but only 
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Msubaru(<l-5 h^o 1 Mpc) [10“ h 70 > M 0 ] 
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Figure 16. Left panel: Black points indicate the mass based on our analysis of CFHT data as a function of the mass obtained from the 
Subaru data. As opposed to the results presented in the right panel, the object catalogs were not matched. We find that the Subaru- 
based masses are 4 =b 6% higher than the CFHT-based result s we report in this paper. The open points indicate the comparison of our 
Subaru-based masses with those from [App legate et al. 1 I20l4 which are 7 it 6% higher. The dotted line indicates the line of equality. 
Right panel: Ratio of the ensemble averaged tangential distortion as a function of radius measured from CFHT and Subaru data. The 
ratio is computed by combining the measurements for the sample of 7 clusters, where the catalogs are matched, such that they contain 
essentially the same objects. The hatched region indicates the 68% confidenc e region for the average ratio (gy FHT /gf, ubaru ) = 0.99±0.02. 
For reference we also indicate the fitting range used in the comparison with lApplegate et al.l (llOlll i. 


present results for an ensemble stacked lensing signal and 
do not provide updated masses for the individual clusters. 

4-1.2 Direct comparison with Subaru data 

The processed Subaru imaging data used by CLASH have 
been publicly releasecjf] for nine of the twenty clusters. We 
retrieved the data for the four clusters that overlap with 
CCCP. To extend the comparison sample, Keiichi Umetsu 
kindly provided us with the data for Abell 209 and Abell 611. 
Observation of Abell 1758 were provided by James Jee. We 
analysed these data using our CCCP weak lensing pipeline. 
We made no modifications and thus assume that our em¬ 
pirical correction for noise bias also applies to the typically 
deeper Subaru data (note that we also apply the size cut 
of rh < 5 pixels in this case). The results presented in Ap¬ 
pendix [D] suggest that our approach, which is a function of 
signal-to-noise ratio and galaxy size, is sufficiently flexible. 

The data that were provided are stacks of dithered ex¬ 
posures. As a consequence multiple chips can contribute to a 
given location, which can lead to a more complex PSF. Our 
observing strategy with Megacam allowed us to avoid this, 
but we note that the same problem occurs for our CFH12k 
data. However, we did not measure noticeable differences in 
the scaling relations based on CFH12k or Megacam data. 
As we did for the CFH12k data, we use the weight images 
to split the data into regions that more or less correspond 

9 http://archive.stsci.edu/prepds/clash/ 


to the chips of the camera, and analyse the resulting images 
using the pipeline described earlier. This analysis is done 
completely independently from the analysis based on the 
CFHT data. Hence we redo the object detection and mask¬ 
ing, identify the stars which are used to model the PSF, 
etc. 

The results are presented in Figure [16] The left panel 
shows a comparison of the weak lensing masses when the 
source catalog is d etermined independent ly from the CFHT 
analysis. Following lApplegate et ahl (12014) . these masses are 
based on the best fit NFW model to the lensing measure¬ 
ments within 0.75 — Sh/^Mpc, adopting a concentration 
C 200 = 4. For most clusters the Subaru data are deeper, 
resulting in a different effective source redshift distribution, 
which we account for. Consequently the correction for con¬ 
tamination by cluster members also differs somewhat. It is 
important to stress this correction works well on average (as 
shown in Figure [T^J, but the statistical uncertainty is larger 
when comparing a small sample of clusters. 

The filled points in the left panel of Figure [16] com¬ 
pare the masses based on the CFHT data (Mnterature) to 
those determined from our analysis of the Subaru data 
(Afsubaru)- We find excellent agreement with a best fit 
Mcfht = (0.97 ± 0.06)ALsubaru* We can use this result as a 
measure of the systematic uncertainty when cluster masses 
are determined independently using different instruments, 
albeit with the same shape measurement pipeline. This com¬ 
parison is made by fitting a parametric model to the lens¬ 
ing signal at relatively small radii, where our correction for 
contamination by cluster members is largest. As discussed 
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in more detail in £14.31 we expect our apert ure masses to be 
more reliable. We find that the masses from I Applegate et al.l 
( 2014 1. indicated by the open points, are somewhat higher. 
We obtain MwtG = (1.07±0.06)Ms u baru. A similar result is 
obtained if we comp are to the six clusters in common with 
llJmetsu et al.l l(2014 ). where we note that our Subaru based 
masses are in excellent agreement (they are ~ 2.4% lower 
on average). 

To examine the performance of the shape measurement 
algorithm further, we created a source catalog where the ob¬ 
jects were matched by position (note that we do apply the 
size cut before matching). Although blending in the inner re¬ 
gions may cause some misidentification, such a comparison 
should eliminate differences in the source redshift distribu¬ 
tion and the contamination by cluster members. The shape 
measurements for individual galaxies are too noisy, and we 
therefore compare the tangential shear profiles from the two 
telescopes. To improve the signal-to-noise ratio even further 
we combine the signals from the seven clusters. This allowed 
us to measure the ratio of the lensing signal as a function of 
distance from the cluster center. The results are presented 
in the right panel of Figure 1161 We find an average ratio 
(l?T FHT /flT Ubaru ) = 0.99 ±0.02, indicated by the hatched re¬ 
gion. This result suggests that our pipeline is able to recover 
the shapes to within 1 ± 2% for different data sets. 

In addition we have matched our CFHT-based measure¬ 
ments to catalogs provided by the WtG team (Von der Lin¬ 
den, private communication). This direct comparison for the 
overlapping sample of clusters showed a remarkable agree¬ 
ment with {gT CGF /gr tG ) = 0.991±0.018. These direct com¬ 
parisons of shear catalogs obtained from observations using 
different telescopes and teams suggest that the pipelines are 
robust. The differences we observe are consistent with the 
statistical uncertainties associated with comparing such a 
small sample of clusters. 


4-1-3 Comparison to Hectospec Cluster Survey 

The infall regions of galaxy clusters provide an interesting 
alternative way to estimate cluster masses at relatively large 
radii. The Hectospec Cluster Survey targeted a sample of 58 
clusters, 14 of which overlap w ith our study. This survey is 
described in lRines et al.l d2013h . who me asured cluster mass 
profiles using the caustic technique fe.g. lDiaferidll999l b We 
compare the estimates f or M 200 from the be st fit NFW mod¬ 
els to those obtained bv lRines et al.l d2013l j. The results are 
presented in Figure 1171 We note that a comparison to the 
velocity dispersions yields similar results. 

The lensing masses are higher than the dynamical 
masses; the average ratio is 1. 22±0.07. If w e ado pt the mass- 
concentration relation from iDuffv et al.l (l2008fl the agree¬ 
ment is worse with an average ratio of 1.33 ± 0.08. Al¬ 
though the comparison sample is small, we observe substan¬ 
tial scatter. The most significant outliers in Figure [T7] are 
MS090 6.5+1110, Abell 1 758 and Abell 2261. In §3.6 of their 
paper, iRines et al. ISqH) comment on indivual clusters, in¬ 
cluding these outliers. The first is part o f a pair of clusters , 
but the infall patterns can be separated. IRines et al.l (l2013l l 
find the higher mass for the other component, Abell 750, 
but comment that MS0906.5+1110 has the higher X-ray lu¬ 
minosity. lOkabe et alj (l2010l i also find a higher lensing mass 
for this cluster, suggesting that the dynamical mass is too 



Figure 17. Comparison of the dynam ical estimates for M 200 
from the Hectospec Cluster Survey from IRines et ahl ll2013l l with 
the best fit NFW model to the lensing data. The dotted line 
indicates the line of equality. 


low. lCoe et al.l d2012i 'l present a detailed study of Abell 2261, 
suggesting that the dynamical mass estimate for Abell 2261 
is low compared to the lensing and X-ray estimates. We note 
the large range in hydrostatic mass estimates for this clus - 
ter, including the measurement from lMahdavi et al.l (|2013l f : 
although lower than the lensing mass , the X-ray mass is still 
larger than the dynamical mass from lRines et al] d2013f ). Fi¬ 
nally Abell 1758 is a mergi ng system, a t the high redshift 
end of the sample studied bv lRines et al .l (12013 1. As a result 
the dynamical mass is not that well constrained. 


4.2 Aperture masses 

Although the NFW profile is a good description for an en¬ 
semble of clusters, as is suggested by stacked weak lensing 
studies (lOkabe et al.l [20131 : llJmetsu et al.ll2014l ). it may not 
be a good model to fit to individual systems. In particular 
the prese nce of substructure m ay lead to incorrect masses 
(see e.g. iHoekstra et al.l 120001 . for a clear example). The 
CCCP sample contains several complex, or merging clus- 
ters, of which Abell 520 has been studied in particular detail 
dMahdavi et al.ir2007i : Ijee et ap2012l . I20l4 l . In these cases a 
more direct estimate of the projected mass, without having 
to rely on a particular profile, would be preferable. 

It is possible to estimate the projected mass within an 
aperture with minimal assumptions about the actual mass 
distribution. However, comparison to other proxies typically 
still relies on deprojected masses, which do depend on the 
assumed density profile. Various esti mators are available , 
but we prefer to use the one proposed bv IClowe et al.l d 1998h : 

[ r 2 0„2 /"'max 

Cc(n) = 2 / dlnr(7t)+ 2 , / d lnr( 7t ), (11) 

Jr 1 'max '2 J r2 

which can be expressed in terms of the mean dimensionless 
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surface density interior to n relative to the mean surface 
density in an annulus from r 2 to r max 

Cc(ri) = R(r' < n) - n(r 2 < r < r max ). (12) 

Hence we can determine the mass up to constant, which is 
determined by the mean convergence in the annulus r 2 < 
r' < r ma x- Assumptions about the mass distribution enter 
in two ways. First of all, we do not measure the tangential 
shear 7 t, but the reduced shear gr = 7t/(1 — «)■ For this 
conversion we use the best-fit NFW model. However, the 
estimate of £ c (r) depends on the lensing measurements at 
large radii and consequently this correction is small. 

The more important dependence on the density profile 
is through the need to estimate the average convergence in 
the annulus. Although the contribution is relatively small if 
we consider large radii, it cannot be ignored for our analy¬ 
sis. H12 used an outer radius r max = 1000” for the CFH12k 
data. However, for these data the azimuthal coverage is in¬ 
complete. We therefore keep the inner radius of 600”, but re¬ 
duce the outer radius to 800”. The annuli are unchanged for 
the Megacam data, i.e. we use r 2 = 900” and r max = 1500”. 

As was done by H12, we use the best-fit NFW model to 
estimate the mean convergence in the annulus. We quantify 
the sensitivity to the profile by varying the normalization of 
the c(M) relation by 20%. This should c apture the variation 
in res ults for different cosmologies (e.g. iBhattacharva et al.l 
l2013l l. The resulting change in projected mass within an 
aperture of l/ijT^Mpc depends on redshift, with a reduction 
(increase) of < 1% for clusters with z > 0.3 if the normal¬ 
ization is increased (decreased). The changes are somewhat 
larger at lower redshifts but at most only ~ 2%. Hence the 
systematic uncertainty in the estimate of the projected mass 
at this radius is remarkably small. 

We do make one small change with respect to 
H12: we include a cont ribution from neighbouring halos. 
lOeuri fe Hamanal (1201 lh show that such a two-halo term 
dominates over the NFW profile on large scales. A conve¬ 
nient way to describe such con tributions is provided by the 
so-called halo model (see e.g., ICoorav fc Shethl 120021 . for a 
review). The implem entation we us e here is described in 
ICacciato et al.1 (l2013h . It was used in ICacciato et al.l (l2014l l 
to model the lensing signal around galaxies. In this context, 
the cluster lensing signal is simply the lensing signal around 
the brightest cluster galaxy of a very massive halo. The con¬ 
tribution from the clustering of halos, the two-lialo term, is 
in its most basic implementation prop ortional to the linear 
matt er power spectrum (but see e.g. Ivan den Bosch et al.l 
l20ll for a more sophisticated implementation). The con¬ 
stant of proportionality is determined by the product of the 
bias of the halo of interest and an average of the halo bias 
over all halo masses in the r ange of interest weig hted by the 
halo mass function (see e.g. ICacciato et aDl2009l l. 

To compute the contribution to the average convergence 
in the annulus we use the functions for the halo mass and 
bias provided by iTinker et al.l (l2010l l. We find that the cor¬ 
rection is small, less than 1% for the radii we are inter¬ 
ested in. This also depends somewhat on the assumed den¬ 
sity profile, but the main source of theoretical uncertainty 
is the halo bias function, especially at the high mass end 
considered here. We gauge the systematic uncertainty in 
the correction by considering diffe rent fitting functions from 
the literature, namely those from ISheth fe Tormenl (11999h ■ 



z 


Figure 18. Ratio of the projected aperture mass within an aper¬ 
ture of radius l/i^Mpc and the mass obtained by H12. The 
hatched region indicates the weighted average ratio (M/Mjji 2 ) = 
1.093 ± 0.016. Note that the error bars indicate only the uncer¬ 
tainty in the updated mass estimates. 

ISheth et al.l d200lf) and ISeliak &; Warrenl d2004l l. Although 
these may differ by up to a factor of a few at the highest 
masses of interest, the corresponding cluster bias varies by 
at most ~ 20% because of the exponential drop-off of the 
halo mass function. We find that the resulting systematic 
uncertainty in the correction itself is ~ 10%, which can be 
safely ignored. 

The resulting projected masses within fixed apertures 
of 0.5 and l/i^Mpc are listed in Table [2] The uncer¬ 
tainties include the contributions from cosmic noise and 
shape noise. Although the variation in the source redshift 
distribution is in pri nciple an additional source of noise 
dHoekstra et al.l 120111 ). the impact is smaller for the larger 
angular s cales used to compute the aperture masses. Fur¬ 
thermore iHoekstra et al.l d201ll ') find that the combination 
of cosmic noise and the variation in the redshift distribution 
leads to a slight reduction in the uncertainty, compared to 
the situation where only cosmic noise is considered. 

In Figure [18] we compare the projected mass to the re¬ 
sults from H12 as a function of cluster redshift. We find no 
significant trend with redshift and obtain a weighted aver¬ 
age ratio (M/A/H12) = 1.093 ± 0.016. The increase in the 
amplitude of the lensing signal because of the correction for 
the effects of noise in the images is partly offset by the in¬ 
crease in the mean source redshift and the change in the 
mass-concentration relation. 

4-2.1 Deprojected masses 

Although the projected masses can be determined robustly 
within an aperture of fixed radius, comparison with other ob¬ 
servations requires the deprojection of the mass estimates. 
To do so we follow the procedure described in iHoekstral 
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Figure 19. Comparison of M|Jq W from the best-fit NFW model 
and the value for M^qq using the deprojected aperture mass. 
The former is based on measurements of the lensing signal at 
radii 0.5 — 2 h^Q Mpc, whereas the latter uses data from radii 
larger than 7*500 which is typically larger than 1 Mpc. As a 
consequence the measurements are almost independent. The red 
dashed line indicates t he line of equa li ty. T he open points indi¬ 
cate clusters for which [Mahdavi et al. I <201311 measured a central 
entropy Kq < 70 keVcm 2 . 

J2007fi . which was also used in H12: at each radius we 
find the NFW model that yields the same projected mass. 
We take the mass of this model, which depends on the 
adopted mass-concentration relation, to be the estimate for 
the deprojected mass. The results for M^ 00 and q 0 , 
using the relation betw een concentration and mass from 
IDutton fe Macci3 (I2014T ). are listed in columns 6 & 8 of Ta¬ 
ble [2] (for reference we also list the correspondin g radius ta). 
We pr esent results for the c(M) relation from Duffy et alj 
(l2008fl in Appendix [Fj 

In figure [TUI we compare M^ 0 to the estimate obtained 
from the best-fit NFW model to the lensing signal at radii 
0.5 — 2/iyQ 1 Mpc. The two measurements are nearly inde¬ 
pendent, because the aperture mass is based on the lensing 
signal at radii larger than r 500 which is larger than lhyg Mpc 
for most of the clusters. From a linear fit to the masses we 
find that /M^ 0 = 0.97±0.03, i.e. the estimates from 

the NFW model are in good agreement on average. 

It is interesting to compare the results for clusters with 
a low and high central entropy Kq, the value of the depro¬ 
jec ted entropy profile at a radius of 2QhjQ kpc, as measured 
bv lMahdavi et akl (l2013ll . The open points in Figure [T9] in¬ 
dicate the clusters with Kq < 70 keVcm 2 , for which we find 
a ratio 1.04 ± 0.06; for the remaining clusters is on 

average 0.95 ±0.04 times smaller than the aperture mass es¬ 
timate. In neither case do we observe a signifant difference 
between the two mass estimates. 

As discussed in H12, the deprojection depends on 
the mass-concentration relation. The results listed in Ta- 



Figure 20. Ratio of the projected aperture mass within an aper¬ 
ture of radius l/i^Mpc when one of the corrections is not in¬ 
cluded, and the final mass estimate. The red histogram shows 
the distribution if we ignore the correction for the multiplicative 
bias in the shape measurement. The black histogram shows the 
decrease in mass if we ignore the contamination by cluster mem¬ 
bers. The blue histogram shows the change in mass if we use the 
source redshift distribution used by H12. 

blc [2] are based on the mass-concentration relation from 
IDutton fc Macc kil d2014j). If we in stead consider the concen¬ 
trations from Duffy et al.l d2008lf . which are ~ 20% lower, 
we find that the resulting M^qq is 8% lower, and Mjq 0 is 
3% lower. Hence, the value for M^ 0 is fairly robust against 
changes in the concentration. 

If we compare to the deprojected masses from H12 we 
find that M^ 0 has increased by 19% on average, whereas the 
increase in M 2500 is 28%. The lar ger change in M 250 0 is du e 
to the higher concentrations from lDutton fe Macciol J2014f) . 
which affect the deprojection. As discussed in H12 (see Fig. 4 
in that paper), lowering the concentration increases the ratio 

mZ w /m^ 0 . 

4.3 Error budget 

Figure l20l shows the distribution of the relative change in the 
projected mass within an aperture of radius l/i^Mpc when 
one of the corrections is not included. The red histogram 
shows that the average mass would be on average reduced 
by a factor 0.84 if we do not include the correction for the 
multiplicative bias in the shape measurement. Ignoring the 
correction for contamination by cluster members reduces the 
mass by a factor 0.90. As indicated by the black histogram, 
this correction varies more from cluster to cluster. This is 
expected, as the correction depends on the cluster redshift, 
the spatial distribution of cluster galaxies and the richness 
of the cluster. If we use the source redshift distribution used 
by H12 the masses would increase by 4%, as indicated by 
the blue histogram. 
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Although each of these corrections is substantial and 
cannot be ignored, they are well-determined. The associ¬ 
ated remaining systematic errors are much smaller than the 
statistical uncertainties for individual clusters. However, the 
precision that is afforded by the full CCCP sample is much 
higher, and we therefore summarize the systematic error 
budget in this section. These are relevant for the discus¬ 
sion in where we examine the scaling relation between 
the lensing mass and the SZ measurements from Planck. 

We assume a 2% systematic uncertainty in the calibra¬ 
tion of the multiplicative bias. Although Figure [9] suggests 
that the empirical correction results in biases < 1% if we 
adopt eo = 0.25, the variation in seeing and the dependence 
of the bias on morphology are expected to lead to addtional 
error. The latter will need to be investigated in more detail, 
but the observed variation in bias as a function of Sersic 
index (see Fig. 0 suggests that the contribution from mor¬ 
phology is ~ 1%. 

As shown by the black points in Figure [12] the mean 
residual contamination by cluster members is at most ~ 2%, 
except for the innermost regions. Although the contami¬ 
nation exceeds this value at some radii for clusters with 
2 > 0.25 (red hatched region), the contamination is still 
smaller than 2% when averaged over the radii of interest. 
For the projected masses, the uncertainty in the extrapo¬ 
lated density profile changes the masses by less than 1%. 

The largest contribution comes from the uncertainty 
in the source redshift distribution, despite the wavelength 
coverage and depth of the COSMOS data: for clusters at 
z = 0.2 the systematic uncertainty in (/?) is ~ 3%, increas¬ 
ing to ~ 8% for the highest redshift clusters in our sam¬ 
ple (2 = 0.55). These estimates include the difference in 
(0 ) between our pho tometric redshift distribution and that 
of lllbert et al.l J2013T ). and the field-to-field variation in the 
four CFHTLS Deep fields. These systematic errors are inde¬ 
pendent from one-another and thus can be added in quadra¬ 
ture, resulting in a total systematic uncertainty of 4.2% at 
2 = 0.2 and 8.5% at 2 = 0.55 for the projected masses in a 
fixed aperture. 

The deprojected masses are more sensitive to the as¬ 
sumed mass-concentration relation. For instance, M^ 0 de¬ 
creases by ~ 3% if we lower the concentr ation by 20% 
(i.e., by switching to the iDuffv et al.l (l2008l l values). This 
shift in concentration is caused by the change in the cos- 
mological parame ters determined by WMAP5, used by 
IDuffv et al.l (1200811. and the more precise Planck values used 
by Dutton fc Mac.c.i^l (120141 1. However, these estimates are 
based on simulations that only include dark matter and thus 
ignore the additional effects of ba ryon phy sics. T he impact 
of this has been studied by |PuffY_et aljj 20ld 'l using hy¬ 
drodynamic simulations. IDuffv et alj ( 201C ) found that the 
change in concentration is < 10% for cluster-mass halos. We 
therefore adopt a similar uncertainty in the concentrations, 
which implies a systematic contribution of ~ 2% for M^ 0 . 

Importantly the systematic errors for the deprojected 
masses are increased because a change in the lensing signal 
also affects the radius corresponding to a particular over¬ 
density. This increases the uncertainty compared to a fixed 
aperture because the enclosed mass increases with radius. As 
a result we estimate a total systematic uncertainty in q 0 
of 6% for clusters at 2 = 0.2, which increases to 12% at 
2 = 0.55. Given the observational cost of calibrating scaling 


relations, it makes sense to avoid introducing such unnec¬ 
essary sources of uncertainty: numerical simulations should 
instead be used to make predictions for the observed lensing 
measurements. 

In 44.1.21 we determined masses for Abell 175 8 and 
the six clusters in common with lUmetsu et aD (l2014l l using 
available Subaru observations. We found that these masses 
were 4 ± 6% higher than our estimates based on the CFHT 
data. Given the differences between the data from the two 
telescopes (e.g., depth, masked areas), the results agree very 
well. Moreover, when we match the object catalogs we find 
that we recover the average tangential distortion within 
1 ± 2%, suggesting our shape measurement pipeline is ro¬ 
bust. 


5 PLANCK SZE SCALING RELATION 


IPlanck Collaboration et alj J2014d l present constraints on 
cosmological parameters using the cluster number counts as 
a function of redshift for a sample of 189 clusters of galax¬ 
ies detected through the SZE by Planck. The estimates for 
the masses from Planck are based on the scaling relation 
between the X-ray hydrostatic mass and Yx, the product 
of the X-ray temperature and the gas mass. This re l ation i s 
calibrated using measurements from lArnaud et al.l (1201014 . 
who studied a sample of 20 nearby relaxed clusters. Conse¬ 
quently a measurement of Yx can be converted into a hydro¬ 
static mass estimate M^jg. These results are then used to 
establish the relation be tween the SZ si g nal F500 and MZv 


(see Appendix A.2.2 of IPlanck Collaboration et al.l 
for details). We denote this hydrodynamic mass estimate as 


Numeri cal _ si mulations (e.g. iRasia et al.l 120061 : 

iNagai et alj 120071 : lLau et aD 120091 ') suggest that such 
mass estimates are biased low. Similarly, iMahdavi et al.l 
(l2013ll studied the scaling relations between X-ray observa¬ 
tions and the weak lensing masses from H12 and found that 
hydrostatic masses underestimate the weak lensing masses 
by 10 — 15% at rsoo- Our updated masses do not change 
this conclusion, and in fact streng then it. In their analysis 
IPlanck Collaboration et al.l (12014a ') assume that the hydro¬ 
static masses are biased low by a factor (1 — b) = 0.8 based 
on a comparison with numerical simulations. They find 
that their best fit parameters for as and are in tension 
with the measure ments obtained from the a n alysis of the 
primary CMB bv IPlanck Collaboration et al.l (l2014bl) . The 
results can be reconciled by considering a low value of 

(1 - b) ~ 0.6. ___ 

Recently, Ivon der Linden et al.l (l2014bll estimated the 
bias using the lensing masses for the 38 clusters in com¬ 
mon between Planck and WtG. They compared their esti- 
mates for M500 b ased on the NFW fits with C200 = 4 from 
lApplega te et al. (|2014|) to the h ydrosta tic mass estimates 
from Planck Collaboration et ahl (l2014al ~). They obtained an 
average ratio (1 — b) = 0.69 ± 0.07, which alleviates the ten¬ 
sion. As our comparison in 44.1.11 and Fig. [TJ] shows, the 
WtG masses are slightly higher than our estimates when we 
follow the same approach, but when we compare the masses 
from WtG to our deprojected aperture masses, which are 
more robust and therefore used here, we find that the agree¬ 
ment is excellent. 
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Figure 21. Left panel : the dep rojected aperture mass M 500 from weak lensing as a function of the hydrostatic mass from 
IPlanck Collaboration et al.| d2014alh Note that MgQQ nck is measured using r $oo from the estimate of Yx, and Mgis determined 
using the lensing derived value for rsoo. The black points show our CCCP measurements, with the filled symbols indicating the clusters 
detected by Planck with a signal-to-noise ratio SNR > 7 and the open points the remainder of the sample. The dashed line shows the 
best-fit power law model. The WtG results are shown as rosy brown colored points. Right panel: ratio of the hydrostatic and the weak 
lensing mass as a function of mass. The dark hatched area indicates the average value of 0.76 ± 0.05 for the CCCP sample, whereas the 
rosy brown colored hatched region is the average for the published WtG measurements, for which we find 0.62 ± 0.04. 


There are 38 clust ers in common between CC CP and 
the catalog provided by lPlanck Collaboration et al.l (|2014al l. 
although we omit Abell 115 from the comparison as we 
determine masses for the two separate components of this 
merging cluster. The left panel in Figure [21] shows the de¬ 
projected aperture mass as a function of the hy¬ 

drostat ic mass M$QQ nck from [Planck Collaboration et alj 
d2014al l. Note that the observed value for Yx was used to es¬ 
timate the radius rsoo used to determine M^ 0 l g nck , whereas 
M500' is based on the v alue for rsoo listed in Table [2] For the 
cosmological analysis, IPlanck Collaboration et~ak d2014d l 
restricted the sample to clusters above a SNR threshold of 7 
in unmasked areas. In our case, the mask only impacts the 
merging cluster Abell 2163, which corresponds to the right¬ 
most point in Figure 1211 There are 20 SNR> 7 clusters in 
common with CCCP and these are indicated as filled points 
in Figure 1211 whereas the remaining clusters are indicated 
by the open points. We find that the SNR threshold is essen¬ 
tially_a 1 _selection_j3y jnass. Fomreference, the measurements 
from Ivon der Linden et all d2014bli are indicated by the rosy 
brown colored points. 

The right panel shows the ratio of the hydrostatic 
masses from Planck and our weak lensing estimates for all 37 
clusters in common. The hatched region indicates our esti¬ 
mate for (1 — 6) = 0.76±0.05 (stat)±0.06 (syst), which was 
obtained from a linear fit to M^ nc k as a function of M§qq 
that accounts for intrinsic scatter ([Hogg et al.ll2Q10l ). The 
systematic error is based on the estimates presented in 11. .31 
We measure an intrinsic scatter of (28±6)%, most of which 
can be attributed to the triaxial nature of dark mat ter ha¬ 
los (e.g. lCorless fe Kfngll2007l : lMenegIietti et al.ll2010ll . If we 


restrict the comparison to the clusters with SNR> 7 (black 
points) we obtain (1 — 6) = 0.78 ± 0.07, whereas (1 — b) = 
0.69 ± 0.05 for the remaining clusters. For reference, the 
rosy brown colored poin ts and hatched region indicat e the 
results for WtG, used in Ivon der Linden et alj (j2014bl l. We 
refit these measurements, which yields (1 — b) = 0.62 ± 0.04 
and an intrinsic scatter of (26 ± 5)%. Our measurement of 
the bias is in agreement with the n ominal value adopted by 
IPlanck Collaboration et alj d2014d l and we conclude that a 
large bias in the hydrostatic mass estimate is unlikely to be 
the explanation of the tension of the cluster counts and the 
pri mary CM B. 

Ivon der Linden et alj d2014bl l find modest evidence for 
a mass dependence of the bias, with Mpi anc k oc AfwtG- It 
is therefore interesting to repeat this for our measurements. 
If we restrict the fit to the clusters with a SNR> 7, the 
range is too small to obtain a useful constraint on the slope. 
We therefore fit a power law to the CCCP measurements of 
the 37 clusters that overlap with lPlanck Collaboration et al.l 
d2014ah . which yields 


YlPlanck 

10 15 /iroM Q 


(0.76 ± 0.04) x 


Mcccp 

10 15 /if o 1 M Q 


0.64±0.17 


and an intrinsic scatter of (21 ± 4)%. The slo pe is similar 
to that found by Ivon der Linden et al.l d2014bh and our re¬ 
sults therefore support their conclusion that the bias in the 
hydrostatic masses used by Planck depends on the cluster 

mass, but our norm alization is 9% h i gher. _ 

As noted above. [Planck Collaboration et all d2014d l use 
X-ray data to relate the observed SZ-signal to cluster mass. 
It is, however, more convenient to directly constrain the 
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scaling relation between the lensing mass and the observed 
SZ signal. H1 2 presented results for the e arly data re¬ 
lease from IPlanck Collaboration et al.l (|201ll ~), but here we 
expand the sample to the 37 clusters in common with 
CCCP and use the measurement s for Y500 provided by 
IPlanck Collaboration et al.l (l2014al ). 

Assuming a constant gas fraction and self-similarity, the 
SZ s ignal Y scales with ma ss as M| qq oc D 2 ng E(z)~ 2 ^ 3 Y 5 oo 
(e.g. iMcCarthv et al.l I 2 OO 3 I ") . where D ang is the angular di¬ 
ameter distance to the cluster and E(z) = H{z)/Hq is the 
normalized Hubble parameter. The results are presented in 
Figure l22l where the open points indicate the clusters that 
Planck detected with a SNR< 7. The dashed line indicates 
the best-fit power law relation to all clusters in common, for 
which we find 


10 15 h7nM E 


0 


(1.01 ±0.06) x 


/ 10 4 D ang Y5oo 
\ v £(2) 2 /3MpC 2 


0.53±0.13 


and an intrinsic scatter of (27 ± 6)%. For this comparison 
we used measurements of the SZ signal and the lensing mass 
in apertures that were determined independently. Although 
this is what one needs for the cosmological interpretation 
of the Planck cluster sample, the use of different apertures 
introduces additional noise as well as an offset because the 
lensing aperture is larger on average. 

In the case of hydrodynamic simulations of clusters the 
comparison can be made at a common radius, as rsoo is 
known. It is therefore useful to consider also the scaling re¬ 
lation for the SZ signal within r^o- The SZ measurements 
within were estimated (Arnaud & Pratt, private com¬ 
munication) using the two-dimensional marginal probability 
distribution between the SZ sig nal and size available from 
the 2 013 Planck SZ catalog (IPlanck Collaboration et al.l 
l2014al ). They correlate very well with Ysoo- The best fit 
power-law scaling relation is given by 


10 15 /2.*7n Mr? 


= (0.98 ±0.05) x 


0 


10 4 Pan g y(r^ O L ) Y 
_E(z) 2 / 3 Mpc 2 J 


and the intrinsic scatter is reduced to (22 ± 7)%. Note that 
for both scaling relations the slopes are consistent with the 
value of 0.6 expected for a self-similar model. 

As a caveat, we note that the combination of relatively 
low significance of the SZ detections and large intrinsic scat¬ 
ter leads to Malmquist bias. As a result the average SZ signal 
of the observed sample is biased high compared to the av¬ 
erage of the parent population the clusters are drawn from. 
Consequently, the best ht parameters for the scaling rela¬ 
tions are expected to be slightly biased. If the CCCP sample 
were a well-d efined sample one could account for this, which 
was done by [P lanck Co llaboration et al.l ( 2014 c). Although 
iMahdavi et al.l (120131 ) show that the CCCP sample is repre¬ 
sentative of more carefully selected samples of X-ray lumi¬ 
nous clusters, the selection may impact the scaling relation. 
We therefore do not attempt to correct our results, but note 
that, b ased on the finding s of Planck Collaboration et all 
(l2014d ) and lMahdavi et al1(l2013l) . we expect the bias to be 
small. 



E z" 2/3D ang Y ( r S)o) [MpC 2 ] 


Figure 22. Plot of M 500 , the aperture mass estimate within an 
aperture , as a function of the SZ signal measured within 
an aperture r^QQ from [Planck Collaboration et all (l2H14al) . The 
dashed line indicates the best fit power law, which has a slope of 
0.53±0.13. We measure an intrinsic scatter of (27±6)%. The open 
circles indicate clusters which Planck detected with a SNR< 7. 


6 CONCLUSIONS 

Accurate cluster masses are necessary to interpret the clus¬ 
ter counts from wide-area surveys. In particular the scaling 
relations and their scatter need to be determined observa- 
tionally. Weak lensing masses are particularly well-suited as 
they provide a direct estimate of the projected mass and do 
not depend on the dynamical state of the cluster. In this 
paper we r evisited the ana l ysis o f a sample of 50 massive 
clusters by iHoekstra et al.l (120121) , with a particular focus 
on improving the corrections for various sources of system¬ 
atic error in the cluster masses. 

We use extensive image simulations to quantify the bias 
in our shape measurement algorithm. The bias is a strong 
function of signal-to-noise ratio and size, but depends rel¬ 
atively weakly on surface brightness profile and ellipticity 
distribution. We demonstrate that the inferred bias depends 
on the input parameters used to create the simulated data. 
For instance we find that the bias is underestimated if faint 
galaxies are lacking from the simulations: the bias converges 
if the simulation includes galaxies that are at least 1.5 mag¬ 
nitude fainter than the limiting magnitude of the sample of 
sources used for the lensing analysis. The large number of 
simulated galaxies enables us to determine an empirical cor¬ 
rection, which is found to be accurate and robust to the main 
uncertainties. We estimate that the systematic uncertainty 
due to the shape measurements is at most 2%. 

The dominant source of systematic error is the source 
redshift distribution, which is needed to convert the lens¬ 
ing signal into an estimate of the mass. We use the 
latest state-of-the-art photometric redshift catalogs that 
are based on measurements in 29 bands in the COS- 
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MOS field dScoville et al.l [20071 : ICapak et a,l.ll2007l "). includ¬ 
ing new deep near-infra red observations from UltraVISTA 
dMcCracken et aLll20120 . Compared to our previous analy¬ 
sis, this leads to a small increase in mean source redshift, or 
a modest (~ 4%) reduction in the cluster mass. Despite the 
unprecedented quality of the data, different analyses show 
variations in the source redshift distributions that result in 
systematic uncertainties that are substantial compared to 
the statistical uncertainties for the full sample of clusters. 

We find that the projected aperture masses within aper¬ 
tures of fixed radius provide the most robust measurements. 
The dependence on the assumed density profile is minimal 
and the systematic error is dominated by the uncertainty 
in the source redshift distribution. We estimate a total sys¬ 
tematic error of 4.2% for clusters at z = 0.2, which increases 
to 8.5% at z = 0.55. To compare to measurements at other 
wavelenghts, we deproject the masses. This leads to an in¬ 
creased sensitivity to the assumed density profile, although 
this is modest in the case of M500: we estimate a systematic 
uncertainty of 2%. Compared to the masses within a fixed 
aperture, the additional uncertainty in determining rsoo in¬ 
creases both the statistical and systematic errors. Compari¬ 
son of the deprojected aperture masses within r 500 and the 
corresponding mass for the best fit NFW profile shows good 
agreement, even though the two estimates are nearly inde¬ 
pendent from each other. 

We compare mass es for the 18 clusters in common with 
lApplegate et al.l d2014l 'l. To do so we fit an NFW model with 
c = 4 to our lensing signal within 0.75 — S/i^Mpc. We find 
that the resulting CCCP masses are on average 8 ± 4 per 
cent lower, which reduc es to 6% if we u s e the same source 
redshift distribution as lApplegate et ahl rt2014h . Given the 
limitations of the comparison, it is not clear whether this 
difference is a sign of residual systematics or merely statis¬ 
tical in nature. 

We did examine the robustness of our analysis. To this 
end we determin ed masses for Abell 17 58 and the six clusters 
in common with lUmetsu et al.l d2014h using available Sub¬ 
aru observations. We found good agreement with our mass 
estimates based on the CFHT data, suggesting our pipeline 
yields robust results when data from different telescopes are 
analysed. Interestingly, a direct comparison of our CFHT 
measurements to catalogs provided by the WtG team (Von 
der Linden, private communication) also showed remarkable 
agreement with (<7 t CCP / gr tG ) = 0.991 ± 0.018. 

The overlap of 37 clusters with the sample of 
clusters for which Planck dete cted the SZ signal 
(IPlanck Collaboration et al.l l2014al l enables us to cal- 
ibrate the bias in the hydro static masses used by 
IPlanck Collaboration et al.l (I2014J 1 to infer cosmological pa¬ 
rameters. The resulting estimates for as and are in 
tension with the measurements f rom the primary CMB 
IPlanck Collaboration et al.l (l2014bl h Our measurements for 
the overlapping clusters yield 1 — b = 0.76 ± 0.05(stat) ± 
0.0 6(syst), in good agreement with th e nominal value used 
by IPlanck Collaboration et al.l (l2014d h Our results do not 
support a large bias in the hydrostatic masses, which could 
alleviate the tension. We also directly calibrate the scaling 
relation between the SZ signal Y500 and the weak lensing 
mass. When we compare the lensing mass to the SZ signal 
measured in the same aperture, we find a best fit slope of 


0.64 ±0.12 which is in good agreement with the ex pectation 
of a self-similar model (e.s. iMcCarthv et al.ll2003l h 

The constraints from the current CCCP sample are al¬ 
ready limited by systematic uncertainties. The most domi¬ 
nant of these is our limited knowledge of the source redshift 
distribution. Although this can be alleviated by measuring 
phot ometric redshifts for the sources in the cluster fields 
(e.g. lApplegate et al .Il201-ll : lUtnelsu et alJl2014l h biases may 
remain due to limited wavelength coverage. Improving this 
situation is critical to calibrate cluster scaling relations to 
the level of accuracy afforded by the next generation of clus¬ 
ter surveys. 
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Figure Al. Top panel: Multiplicative bias as a function of the 
ellipticity of the PSF for galaxies with 22 < m r < 25. The red 
hatched region indicates the results obtained for a circular PSF. 
The histogram indicates the distribution of PSF ellipticities ob¬ 
served in the CFHT data used in this paper. Bottom panel: The 
average shear 71 measured after PSF correction as a function 
of PSF ellipticity. The average should vanish if the correction is 
perfect. These results suggest that the residual is ~ 4% of the 
original PSF ellipticity. 
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APPENDIX A: EFFECTS OF PSF 
ANISOTROPY 

To quantify the multiplicative bias in the shear measure¬ 
ment resulting from noise in the images and limitations in 
the correction for the size of the PSF and the weight func¬ 
tion, we created simulated images with a circular PSF. In 
real data, however, the PSF is generally anisotropic, albeit 
to varying degree. PSF anisotropy leads to additive biases 
by introducing coherent alignments in the observed shapes. 
In the case of cosmic shear this is a dominant source of sys¬ 
tematic, but in our case the signal is averaged azimuthally 
around the clusters, which also averages out most of the PSF 
anisotrop\P1. 

Nonetheless it is important to examine whether the re¬ 
sults obtained using circular PSFs can be applied to our 
CFHT data. We therefore created a set of simulations in 

10 We note that this is not always the case, as some imagers show 
strong radial patterns, which can bias the azimuthally averaged 
lensing signal if the cluster is located at the center of the field-of- 
view. 
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Figure Bl. Multiplicative bias for galaxies with 22 < m r < 
25 as a function of the density of stars with 18 < m r < 20, 
based o n the Besanqon mod el of stellar population synthesis of the 
Galaxy <1 Robin et al.|[2003h . The red hatched region indicates the 
results for the simulations without stars. The arrows indicate the 
corresponding Galactic latitudes (for the longitude of Abell 1835, 
l = 340°). The histogram shows the distribution of star densities 
of the CCCP data. 

which the PSF is elliptical (along the x— axis) using our fidu¬ 
cial ellipticity distribution (eo = 0.25). These images were 
analysed as usual, and the resulting multiplicative bias as a 
function of PSF ellipticity is presented in the top panel of 
Figure IaTI Even for rather cllongated PSFs (e > 0.1) the 
increase in bias is at the percent level at most. For compar¬ 
ison, the histogram shows the observed distribution of PSF 
ellipticity in the CCCP data. These results suggest that we 
do not need to account for the PSF anisotropy explicitly in 
order to quantify the multiplicative bias in our data. 

We do not expect the correction to work perfectly, ei¬ 
ther due to limitation of the correction method itself, or due 
to the fact that the noise in the images biases the polariz¬ 
abilities. The bottom panel shows how well the correction 
for PSF anisotropy performs. We observe a linear trend of 
the recovered average shear as a function of PSF ellipticity. 
For the galaxies with 22 < m r < 25 this results in a residual 
bias of ~ 4% of the original PSF ellipticity. Although the 
bias is smaller for bright galaxies it does not vanish, sug¬ 
gesting that a large part of the observed bias is due to a 
fundamental limitation of the KSB methoTI. 


APPENDIX B: CONTAMINATION BY STARS 

The simulated data that were used to study the multiplica¬ 
tive bias contained a small number of bright stars, which 

11 The KSB algorithm assumes that the PSF is described as the 
convolution of a circular kernel with a compact anisotropic one. 
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were only included to keep track of the PSF properties. In 
real data, however, stars may blend with galaxies or other 
stars. As a consequence they might be misidentified and be 
included in the galaxy catalog, contributing to multiplica¬ 
tive bias (as they are not sheared). 

In this section we study how the multiplicative bias de¬ 
pends on star density by including stars in the images with 
a realistic number density and magnitude distribution. We 
use the Besangon model of stellar populat ion synthesis of 
the Galaxy for the stars (iRobin et al.ll200 3 : ). As a reference 
we consider again Abell 1835, with (1,6) = (340,60), where 
we change the Galactic latitude b to change the star number 
density. We place the stars at random locations in the im¬ 
ages, with to = 0.25, and proceed as before with the shape 
analysis. 

The results for galaxies with 22 < m r < 25 are pre¬ 
sented in Figure iBll which shows the multiplicative bias as 
a function of the density of stars with 18 < m r < 20. For 
densities n s tar < 1.5 arcmin -2 the bias is consistent with the 
results without stars (indicated by the red hatched region). 
For the adopted longitude of l — 340° this corresponds to 
b > 35°. For lower Galactic latitudes, the higher number 
density of stars can lead to appreciable levels of bias. The 
histogram in Figure [Bll shows the distribution of star den¬ 
sity in the CCCP data, which suggests that we can safely 
ignore the contribution from stars. 


APPENDIX C: DETERMINING THE 
EMPIRICAL CORRECTION 


In this appendix we describe the empirical correction used in 
the analysis of the CCCP data. We assu me that the bias is 
a fun c tion of the signa l -to-no ise ratio (e.g. iMelchior fe Violal 
120121 ; iKacprzak et alJ 1201 2) and depends on how we 1 a 
galax y is resolved (e.g. Massev et al.l 120131 : iMiller et al.l 
l2013h . That these parameters are important is also sug¬ 
gested by Figures [o] and [ 6 j which show a larger bias for 
fainter galaxies and larger PSFs. 

Our implementation of KSB provides an estimate for 
q e , t he uncertainty in the polarization (see iHoekstra et al.l 
l2OO0l. for details). The reciprocal of this quantity is a useful 
proxy for the signal-to-noise ratio v, hence we use v = 1 /<r e . 
To quantify how well a galaxy is resolved we use 


iz 2 = 


rl, 


.2 _ 2 ’ 

h , gal ' h, 


(Cl) 


where rh,* denotes the half-light radius of the PSF and rh , ga i 
that of the observed galaxy. The denominator corresponds to 
the unconvolved size of a source if galaxies were Gaussians. 
Importantly, these are quantities that can be measured for 
individual sources. 

We consider values of eo between 0 and 0.5 with steps 
of 0.05 and create 169 pairs of images, each 10,000 by 10,000 
pixels, with constant shears of —0.06 < 7 i < 0.06 applied. 
As a result, for each eo we analyse ~ 10 7 galaxies with 20 < 
m r < 25 (note that the input catalog does contain fainter 
galaxies). For a given eo we bin the measurements in fine 
bins of v and rh and determine the bias. 

We considered various fitting functions, with the aim to 
find an adequate correction that is also robust against the 


uncertainty in the true value of eo- We consider a conserva¬ 
tive range of eo 6 [0.15, 0.3], where we note that a stronger 
prior on the input ellipticity distribution would allow the 
uncertainties to be reduced, and the fitting functions to be 
optimized. 

Although further optimization is possible by introduc¬ 
ing additional parameters, we opted for a correction with 4 
free parameters: one to describe the size dependence of the 
bias for a given v, and three to capture the dependence on 
v. The correction takes the form 


1 + TV) 


b(v) 

1 + a(eo)7 Z' 


(C2) 


The first step is to determine a(eo), which is done by exam¬ 
ining fj, as a function of size for narrow bins in v. The left 
panel of Figure [Cl] shows the measurements as a function of 
half-light radius for galaxies with 20 < v < 30 for three ellip¬ 
ticity distributions. The number of objects as a function of 
observed half-light radii is shown in the bottom panel; most 
sources are small (the value for the PSF is rh,* = 2.056 
pixels). 

For each ellipticity distribution we determine the best 
fit value for a, under the assumption it only depends on eo- 
The results are presented in the right panel of Figure [ClI and 
listed in Table ED The right panel in Figure ED shows that 
a(eo) increases smoothly with increasing eo. The lines in the 
left panel indicate the predicted bias, with the amplitude a 
free parameter (which is used to determine b(u ), see below). 
Our parametrization for the size dependence does fairly well 
for the bulk of the sources, but it does not capture the results 
for rh > 4 pixels. Better agreement is obtained if we include 
an additional term oc TZ 2 , but we found that this did not 
improve the robustness of the correction. Similarly, we found 
that we could have included a dependence on 1/v, again with 
limited effect. Note there is covariance between some of the 
parameters. For instance, for small values for v we do expect 
b(y) and a(eo) to be correlated because those galaxies have 
on average larger values for TZ. 

Closer inspection of the objects with large observed 
sizes revealed that most of these are faint objects for which 
the input sizes were much smaller or blends with other galax¬ 
ies. As the latter might be particularly relevant for the study 
of galaxy clusters, we decided to only use galaxies with 
rh < 5 pixels for the actual lensing analysis. The ensemble 
averaged lensing signal as a function of source size presented 
in Figure [TT] indicates that the results for large sources are 
indeed biased low. Note that we do not apply this cut when 
fitting for a(eo), because we found that the correction per¬ 
formed a bit better when we considered the full range in 
sizes. 

The next step is to quantify how fi depends on the 
signal-to-noise ratio. The results for eo = 0,0.25 and 0.5 
are presented in the left panel of Figure ED For low eo the 
bias increases monotonically, asymptoting to a value ~ 0.95 
for large u, whereas for eo = 0.5, b(v) increases first before 
declining. To capture the variation in b(y) we adopt 


6 M = 6o(eo) + ^ + ^). (C3) 

We note that fixing bo ~ 0.95 also gave reasonable results. 
Furthermore 61 and 62 are highly (anti-)correlated, and it 
thus might be possible to reduce the number of free pa- 
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Figure Cl. Left panel: bias 1+^asa function of half-light radius for galaxies with 20 < v < 30. The black points indicate the results 
for our adopted value of eo — 0.25; the red (blue) points are for the extreme values of eo = 0 (eo — 0.5). The lines indicate the fits using 
our best fit value for a(eo). The bottom panel shows the distribution of half-light radii for these sources. Right panel: Resulting best-fit 
value for a(eo) for our simulated CCCP data. The values are listed in Table eh 
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Figure C2. Left panel: b(v) for the three ellipticity distributions, as well as the best fit models. The bottom panel shows the distribution 
of signal-to-noise ratios for galaxies with 20 < m r < 25. Right panel: the model parameters bi as a function of eo, which show a smooth 
dependence. The values are listed in Table [Cl] Note that bo is fairly close to unity. 


rameters in principle. Interestingly the value for bo, which 
corresponds to the bias for bright, large sources, is fairly 
close to unity. We list the best fit parameters in Table ED 
and the right panel in Figure EH shows that the parameters 
vary smoothly with eo- 


In i)2.4l we discuss how well the correction performs as 
a function of ellipticity distribution and seeing. We also ex¬ 
amined Hcoi as a function of apparent magnitude and found 
that residuals are < 1% for galaxies fainter than m r = 22, 
whereas /r c or ~ 0.02 for galaxies with 20 < m r < 22. 
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Table Cl. Parameters for the empirical correction for multiplicative bias 


€ 0 

a 

bo 

bi 

&2 

0 

0.056 ±0.0015 

0.952 ± 0.0143 

0.84 ±0.15 

-5.50 ±0.32 

0.05 

0.067 ±0.0017 

0.957 ± 0.0145 

0.86 ±0.15 

-5.47 ±0.33 

0.10 

0.107 ±0.0022 

0.958 ± 0.0164 

1.20 ±0.17 

-6.15 ±0.37 

0.15 

0.164 ±0.0029 

0.959 ± 0.0233 

1.61 ±0.24 

-6.81 ±0.53 

0.20 

0.247 ±0.0038 

0.970 ± 0.0235 

2.01 ±0.24 

-7.34 ±0.53 

0.25 

0.348 ± 0.0051 

0.988 ±0.0217 

2.29 ±0.22 

-7.40 ± 0.49 

0.30 

0.473 ± 0.0067 

0.993 ± 0.0220 

2.73 ±0.22 

-7.72 ±0.50 

0.35 

0.616 ±0.0087 

1.017 ±0.0298 

3.04 ±0.30 

-7.70 ±0.68 

0.40 

0.729 ±0.0104 

0.997 ± 0.0260 

3.57 ±0.26 

-8.32 ±0.59 

0.45 

0.864 ±0.0120 

1.018 ±0.0289 

3.85 ±0.29 

-8.43 ± 0.65 

0.50 

0.921 ±0.0121 

1.002 ±0.0265 

3.95 ±0.27 

-8.18 ±0.60 


This overcorrection is probably caused by the fact that our 
parametrization tends to underestimate the bias for large 
objects. 

The bias depends on Sersic index (see Fig.[H|), and there¬ 
fore the performance of the empirical correction may differ. 
Figure [C3l shows the residual bias as a function of elliptic- 
ity distribution for four different Sersic indices. The black 
points with errorbars indicate the results for the ensemble 
of galaxies, with half-light radii and Sersic indices drawn 
from the observed distribution of values. The lines indicate 
the results for the simulations where the Sersic indices were 
fixed to the indicated values. The range in bias is similar to 
what was observed in Figure [HJ with a positive bias in the 
expected range of eo for n = 4 and a negative bias for n = 1. 
Given the relatively weak dependence of Sersic index and the 
fact that the distributio n of Sersic values is well determined 
from HST observations (iRix et al.l l2004| j , we conclude that 
our correction for the ensemble is robust. Hence, we ignore 
the dependence of the bias on Sersic index. 


APPENDIX D: APPLICATION TO SIMULATED 
RCS2 DATA 

Our KSB implemen tation has also been used by 
Ivan IJitert et al J (1201 il l to measure the lensing signal around 
galaxies using data from the second Red-sequence Cluster 
Survey (RCS2). It is therefore interesting to examine the im¬ 
pact of our findings on those results. Compared to our CCCP 
data (with a total integration time of 3600s per position) the 
RCS2 data are much shallower, consisting of a single 480s 
exposure in the riband. We create a separate set of simula¬ 
tions where the noise level matches that of the RCS2 data. 
The res ulting bias for galaxies w ith 22 < m r < 24 (the range 
used by Ivan Uitert et al.l d201lh for the RCS2 analysis) as a 
function of eo is presented in the lo wer panel of Figu re ID II 
The actual bias is smaller because Ivan Uitert et al] (120110 
used the STEP1 implementation, also used in the CCCP 
analysis of H12. We also note that Ivan Uitert et ah Mm3) 
used SExtractor to detect objects, and used the resulting 
half-light radius for the weight function in the subsequent 
shape analysis. 

Here we are interested whether our correction scheme 
can be used to data with different noise properties. We there¬ 
fore apply the correction to the simulated RCS2 results and 
present the residual bias in the top panel of Figure [Dl] The 



Figure C3. Bottom panel: Residual multiplicative bias for galax¬ 
ies with different values for the Sersic index as a function of eo- 
The simulated galaxies have the same distribution in half-light 
radii as the regular simulations (indicated by the black points 
with errorbars) and 22 < m r < 25. As in the real data we only 
include galaxies with an observed size < 5 pixels. The full 
range spans about 0.03, but the uncertainty for the ensemble of 
sources is much less as the distribution of Sersic profiles has been 
obtained from HST observations. 

black line indicates /i cor if we use our reference value of 
eo = 0.25 to correct the measurements for the various ellip- 
ticity distributions. The results suggest that the parameters 
that were optimized to correct CCCP can be used for the 
shallower RCS2 data as well. 


APPENDIX E: OBSCURATION BY CLUSTER 
MEMBERS 

As described in mm we account for the dillution of the lens¬ 
ing signal due to cluster galaxies in the source galaxy sample 
using the excess galaxy counts as a function of cluster-centric 
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Figure Dl. Bottom panel: Multiplicative bias for galaxies with 
20 < m r < 24 as a function of eo for a simulation of RCS2 data, 
which are shallower than the CCCP data. Top panel: Residual 
bias after we apply our correction, using parameters optimized 
for CCCP data, to the simulated RCS2 data. The solid line indi¬ 
cates ficor if we use eo = 0.25 to correct the simulations for other 
ellipticity distributions, suggesting our approach is adequate for 
these data as well. 



r / r 500 


Figure El. Plot of the correction of the counts of galaxies with 
22 < m r < 25 for the obscuration by cluster galaxies as a function 
of radius in units of rsoo- The black points indicate the average 
from a representative subset of clusters, whereas the light grey 
points show the individual measurements. The red line is the best 
fit model as described in the text. 


radius. Figure [TT] shows that the number density of bright 
cluster galaxies is substantial at small radii, which may affect 
t he de tection of sources. As shown in lSimet fe Mandelbauml 
d2014h . this is an important source of bias for the measure¬ 
ment of the magnification signal, but also is relevant here 
as it leads to biases in the correction for contamination by 
cluster galaxies. 

To quantify the impact of the obscuration by cluster 
galaxies on the source galaxy counts, we use our image simu¬ 
lations: we add the cluster observations to the simulated im¬ 
age and perform the object detection and analysis on the im¬ 
ages with and without the cluster addecP^I. We identify the 
objects detected in both catalogs and measure their number 
density as a function of cluster-centric radius. As expected, 
at small radii we observe a decrease in the recovered number 
density. 

Figure EO shows the corresponding correction for the 
source counts as a function of radius in units of rsoo. The 
black points correspond to the average from a representative 
sample of clusters that we used in this study. The individ¬ 
ual measurements are indicated by the lightgrey points. For 
radii larger than rsoo the observed excess counts are biased 
low by a few percent. Hence, our aperture mass estimates for 
M 500 are unaffected. On the smallest scales considered for 
the NFW fits, the observed counts are biased low by ~ 10%. 
Even in this case the impact is small, as this is a 10% correc¬ 
tion to a correction that itself is 30%; the resulting change 
in the best fit mass is a ~ 1 — 2% increase. We find that the 
correction can be described by 

_L_ 0.021 

/obscured 0.14 + (r /r$00 ) 2 ’ 

which is the red line shown in Figure EO We use this model 
to correct the measurements of the contamination by cluster 
members. 


APPENDIX F: MASS ESTIMATES USING 
DUFFY ET AL. (2008) C(M)-RELATION 


Previous cluster weak lensing studies, including H12, pre¬ 
sented mass measureme nts using the relati on between mass 
and concentration from iDuffv et al.l (120081 1 , which is based 
on n umerical simulatons assuming a WMAP5 cosmol¬ 
ogy dKomats u et al.l 2009ll. The first res u lts from Planck 
presented by [Planck Collaboration et ahl d2014bl l suggest 
higher values for both the normalization of the matter power 
spectrum as and the mea n density ff m . For t his re ason we 
adopted the relation from iDutton fc Maccibl d2014|) , which 
yields a concentration that is ~ 20% higher than Duffy et al.l 
(1200811 for a given mass. 

This affects the masses inferred from parametric NFW 
fits and the deprojected aperture masses, but not neces¬ 
sarily in the same sense. For instance, if we switch to the 
IDuffv et al.l (|2008l ) relation, the estimate for M 2500 decreases 
by ~ 7%, no matter whether we consider the aperture mass 


12 We ignore magnification which increases the fluxes of sources 
and thus reduces the effects of obscuration somewhat. Note that 
the number density of sources is not affected significantly by 
magnification because the power law slope of the number counts 
dlog 7Vg a i/dM ~ 0.38 — 0.4, as discussed in J3] 
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Table FI. Weak lensing mass estimates for the CCCP sample using the mass-concentraton from iDuffv ct al.l d'2CjO>-j ) 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

name 

MET 

MfT 

712500 

1+500 

7*500 

MlEo 

m nfw 

vir 

m nfw 

iW 2500 

M 5 Nfw 

Abell 68 

5.0 ±0.4 

8.7 ± 1.2 

560 

q O+0.3 
-0.3 

1391 

9-9+6 

14 2"*” 3 ' 1 

1 ^' z — 3.1 

9 q+0.6 
^•°-0.6 

7-8+t 

Abell 209 

3.7 ±0.5 

7.8 ± 1.5 

508 

9 q+0.4 
Z ‘°- 0 A 

1249 

6-8±i;« 

9-31+ 

1 9+0-5 
1,y -0.5 

5 1+ 1 - 5 

Abell 267 

4.6 ±0.5 

7.8 ± 1.5 

546 

2 9 +0 ' 4 

z ‘ y -0.4 

1222 

6-61+ 

8.6+8 

1 fi+0-6 

^•0-1.4 

Abell 370 

8.1 ±0.6 

16.9 ± 1.7 

661 

6 1+ 0 - 6 
O1 -0.6 

1659 

19-3+5 

34.61®;® 

6.2I+ 

18.7+3 4 

Abell 383 

2.8 ±0.6 

6.4 ± 1.4 

445 

1 3+0.5 

1232 

6 4 +2,2 
u -2.4 

6-2+6 

1 q+0.6 
1,O -0.6 

3-51+ 

Abell 963 

4.5 ±0.5 

7.9 ± 1.4 

555 

o q+ 0.5 
0 - u -o.5 

1274 

7.21+ 

13.6+;® 

9 7+0.7 

Z '' —0.7 

7-5+1 

Abell 1689 

8.6 ±0.5 

15.0 ± 1.4 

734 

5 7+°- 6 
°'-0.6 

1616 

14 -4 +1 

35.0+° 


18.51+ 

Abell 1763 

5.4 ±0.5 

10.5 ± 1.4 

615 

4 1+0-6 
* -0.5 

1529 

i 9 7+3.3 
iZ *'-2.9 

19.0+1 

0 0+0.8 

^-o.s 

10.3+1 

Abell 2218 

5.7 ±0.5 

9.3 ± 1.4 

646 

4 6 +0 ' 7 
^• d -0.6 

1403 

9 4+ 2 - 2 
y, ^-2.2 

18.31+ 

q c+0.8 
o,o -0.8 

9 9 -1 " 2 ' 4 
y,y -2.3 

Abell 2219 

4.5 ±0.6 

10.1 ± 1.4 

562 

q 9 + 0.6 
°‘ z -0.6 

1418 

10.21+ 

12-5+1 

9 c+0.6 
z,o -0.6 

6-9+1 

Abell 2390 

6.1 ±0.5 

11.6 ± 1.3 

648 

4 g+0-6 

0.6 

1407 

10-0+5 

26.4+;® 

4 9+0- 8 
^-0.8 

14-2+3 

MS 0015.9+1609 

8.5 ±0.6 

20.8 ±2.0 

644 

6 9 +0 ' 7 

1654 

23-41 ®;£ 

33 4+ 8 - 5 
00, ^_ 8.0 

o.y_i 4 

18.41+ 

MS 0906.5+1110 

4.5 ± 0.5 

9.6 ± 1.3 

574 

q 9 + 0.5 
-0.5 

1457 

10-41+ 

14-2+1 

9 0+0.6 
^•°-0.6 

7 7+1-8 
*•* -1.8 

MS 1224.7+2007 

2.0 ±0.6 

2.6 ± 1.6 

359 


794 

2 0 +0 - 8 
Z,U -0.7 

5-2+1 

1 ^0.6 

1 -0.5 

3 0”*” 1 ' 5 
J,u -1.3 

MS 1231.3+1542 

1.2 ±0.4 

0.5 ± 1.2 

359 

0-8+1 

584 

0 7 +oa 

u - 1 -0.4 

1-9+1 

Q c+0.3 
U '°-0.3 

1 1+O.8 

1 ' 1 — 0.7 

MS 1358.4+6245 

4.8 ±0.7 

9.4 ± 1.7 

542 

q 9+0.7 
°‘ Z -0.7 

1316 

9 i+ 2 - 2 
-2.2 

15-0+1 

2 9 +0 ' 7 
z,y -0.7 

Q q+2.1 
° ,ci -2.0 

MS 1455.0+2232 

4.3 ±0.4 

8.1 ± 1.3 

538 

9 g+0.4 
Z ' y -0.4 

1191 

6 2”*” 1 ' 3 
D.Z-1 2 

14 6“*" 2 ' 8 
h,u _2.8 

2 9+ 0,6 
z,y —0.6 

8-0+6 

MS 1512.4+3647 

1.5 ±0.6 

4.5 ± 1.6 

274 

n 4+0- 3 

U -^-0.3 

829 

2 4+ 1,8 

4.1+;? 

0 9” 1 ” 0,4 

u,y -0.4 

2-41+ 

MS 1621.5+2640 

5.4 ±0.6 

12.3 ± 1.6 

565 

4 O" 1 ” 0,7 

^• U -0.7 

1301 

9 Q+ 20 
y - y -1.8 

21.8+1 

4 o +0 - 9 
^• U -0.9 

12-0+1 

CL 0024.0+1652 

6.7 ±0.6 

12.5 ± 1.8 

597 

4 6 +0 ' 6 
^°-0.6 

1357 

10 7 +1 - 9 
1U *‘ -1.9 

97 7+6.8 

Z ‘ * —6.4 

3 i+l- 2 
°* 1 —1.2 

15-1+6 


Abell 115N 

Abell 115S 

Abell 222 

Abell 223N 

Abell 223S 

Abell 520 

Abell 521 

Abell 586 

Abell 611 

Abell 697 

Abell 851 

Abell 959 

Abell 1234 

Abell 1246 

Abell 1758 

Abell 1835 

Abell 1914 

Abell 1942 

Abell 2104 

Abell 2111 

Abell 2163 

Abell 2204 

Abell 2259 

Abell 2261 

Abell 2537 

MS 0440.5+0204 

MS 0451.6-0305 

MS 1008.1-1224 

RX J1347.5-1145 

RX J1524.6+0957 

MACS J0717.5+3745 

MACS J0913.7+4056 

CIZA J1938+54 

3C295 


1.4 ±0.5 

5.3 ± 1.3 

283 

0 4+0.3 
u ’^-0.4 

1087 

4 4+ 1 - 3 

^—1.9 

6-6+1 

1 4+0-5 
1 -^-0.5 

q 7+1-2 

2.6 ±0.5 

6.0 ± 1.3 

399 

l 1+0-5 
i ‘ ± -0.5 

1116 

4-8+3 

7-9+1 

l 0+0-5 

4 4+ 1 - 3 
^-1.3 

2.9 ±0.4 

6.9 ± 1.2 

433 

1 4+0-5 

1 -0.7 

1165 

5-6+;? 

7-1+1 

1 3+0.5 
4 '0 —0.5 

4 0 +1 ' 4 

3.0 ±0.5 

7.6 ± 1.2 

448 

1 0+0-5 
i -°-0.7 

1226 

6.5+3 

10.0+° 

2 O” 1 ” 0-6 
z,u —0.6 

3 3+1-7 
°-°-1.6 

3.0 ±0.4 

8.4 ± 1.2 

451 

1 0+0-5 

1355 

q 7+I.6 

q q+2.9 

5.0_2 7 

1 q+0-6 
1 -®-0.6 

4 9 +1 ' 6 

3.6 ± 0.5 

7.4 ± 1.4 

513 

9 q + 0.5 
Z,O -0.5 

1201 

6 0 +1,5 
uu —1.3 

16.9+-® 

q q+0.7 
O,O -0.7 

9 2”*” 1 ' 9 
y ‘ z —1.8 

3.1 ±0.5 

9.0 ± 1.4 

401 

1 2 +1 ‘ 0 
1 - Z - 0.5 

1321 

8-5+1 

1 9 9+3.5 
1Z,Z —3.4 

2 4+0- 7 
Z -^-0.7 

6-71+ 

2.5 ±0.6 

6.1 ± 1.4 

430 

1 q+0.5 

1203 

5-9+2 

5-0+1 

: 1 +0-5 
l - L - 0.5 

2.8+1 

3.7 ±0.5 

9.0 ± 1.3 

489 

9 9+O.6 
Z -0.5 

1226 

7-1+1 

101+ 3 - 4 

1U - 1 -3.1 

2 0 +0 ' 7 
z,u —0.6 

5.6+1 

4.6 ±0.5 

10.5 ± 1.4 

551 

q 9 + 0.5 
°‘ z —0.5 

1417 

io-8±|;S 

1 3 7+3.8 
4°./_3 6 

0 1+0.7 
°" L —0.7 

8.6+1 

5.4 ±0.5 

12.2 ± 1.4 

540 

0 4+0.4 
^‘^—0.4 

1348 

10.7+j'® 

24 l+ 6-2 
z ^' 1 -5.7 

4 4+ 1 - 1 

^-1.1 

13-2+1 

5.0 ± 0.5 

10.9 ± 1.4 

580 

q 7+O.6 
° m ‘ -0.6 

1333 

9 0 +1,3 

y - u -1.5 

22.3+I 

a 9+0- 8 
^■ z —0.8 

12 1+ 2 - 4 
1Z,i —2.4 

2.5 ±0.5 

4.4 ± 1.4 

439 

1 4+0- 3 

982 

3.2+1 

8-2+1 

1 7+0.5 
—0.5 

45+ 1 . 4 

2.7 ±0.5 

5.7 ± 1.1 

426 

! 3+0-4 
1,Ci -0.5 

1082 

4-4±8;® 

9-71+ 

2 O -1-0 ' 5 

z - u —0.5 

5 4+ 1,5 

5.5 ±0.5 

12.2 ± 1.4 

633 

4 g+0-6 
^ -0.6 

1491 

12.6+1° 

21 6" 1 ” 3 ’ 6 
zl ‘°-3.6 

4 l +0 - 7 
^• 1 -0.7 

11 7+ 20 
±1 - ‘-2.0 

5.3 ±0.4 

10.7 ± 1.2 

603 

4 0 + °- 4 
^• u -0.4 

1387 

9-8+s 

22.4+I 

A 9+0- 8 
^' z —0.8 

12.1+I 

3.7 ±0.5 

7.9 ± 1.2 

522 

2 4+ 0-4 
z ‘^-0.4 

1287 

7 2+ 1 - 4 

' * z —1.3 

15.1+1 

O 0+0.5 
^^-0.5 

8.2+I 

3.8 ±0.5 

7.6 ± 1.2 

519 

9 c+0.4 
z,o -0.4 

1204 

6 2” 1 ” 1 ' 3 
°- z -1.2 

15-1+1 

0 0+0.6 
°- u -0.6 

8.2I+ 

4.2 ± 0.5 

10.3 ± 1.2 

583 

q q+0.6 
0, °—0.6 

1422 

9 5 +1 ' 9 
y, °-1.6 

17-8+1 

0 4+O.8 
o,f *-0.7 

9-6+1 

3.9 ± 0.5 

6.6 ± 1.3 

517 

9 c+0.4 
z,o -0.4 

1156 

3 3+1.6 

10.1+! 

9 1+0.5 
Zll -0.5 

5-6+1 

4.4 ±0.5 

9.5 ± 1.4 

562 

q 1 +0.6 

•^-o.s 

1456 

10.8++ 

19.6+1 

q 7+0.8 
' —0.8 

10.6+1 

4.8 ±0.5 

11.0 ± 1.0 

619 

0 q+ 0.5 
°- y -0.5 

1490 

10 - 9+1 

22 4“*" 3 ’ 8 

4 2” 1 ” 0-7 

^' Z —0.7 

12.0+g 

2.4 ±0.5 

5.6 ± 1.3 

417 

1 2 +0 - 4 
1 ‘ z —0.4 

1106 

4-5+1 

q q+3.1 

».0_2 7 

1 q+0.6 
—0.6 

4-8+1 

6.0 ±0.4 

14.2 ± 1.3 

666 

r q+0.5 
°- c> -0.5 

1654 

16 - 1+1 

28.1+1 

c 9+O.9 
OlZ -0.9 

15-1+1 

5.4 ±0.6 

10.1 ± 1.4 

585 

3 g+0.6 
-0.6 

1302 

8-5+1 

23.6+1 

4-41+ 

12-9+1 

2.9 ±0.5 

2.8 ± 1.3 

457 

1 6 +0 - 4 
i -°-0.4 

890 

2 4+ 0 - 7 
Z ‘^—0.7 

q 3+2.O 
0 - 0 -1.8 

0-81+ 

2-0+1 

4.4 ±0.7 

8.4 ±2.1 

453 

9 4+O.6 
Z ^-0.6 

1071 

6 3 +1 ' 7 

°-°-1.9 

19-4+1 

3.6+1 

10-8+1 

4.1 ±0.4 

8.2 ± 1.4 

507 

9 c+0.4 
z>o -0.3 

1168 

6 2”*” 1 ' 3 
°- z -1.2 

18.2+1 

q c+0.8 
O.0_o 7 

10.0+1 

5.2 ±0.8 

10.1 ± 1.9 

530 

0 4+O.9 
°-^-0.8 

1301 

10 i+ 3 - 2 

iU-1 -3.0 

22.0+1 

4 1+ 11 

12.2+1 

2.1 ±0.9 

6.5 ±2.1 

245 

q 4+O.4 
u ‘^—1003.0 

961 

4 4+ 2 - 2 

^•^-1.9 

7 1+4.8 

' - — 4.2 

1 4+ 1 - 0 
X - 4 -0.9 

4-1+1 

6.4 ±0.9 

19.3 ±2.3 

586 

5 2 4-1 ' 1 
°- z -3.1 

1470 

16 4 -1-3 ’ 4 

44 8+ 10 - 8 
^•o_10.3 

7 6"*" 1,8 
' 1.8 

24.4+1 

3.1 ±0.5 

5.3 ± 1.5 

380 

1 9+O.5 
■ L ‘ Z —0.5 

935 

q 7+1-0 

1 -0.9 

6-9+1 

1 4+0-8 
1, ^-0.7 

4 - 0+1 

5.3 ±0.5 

11.3 ± 1.3 

588 

q q + 0.4 
^•°-0.4 

1557 

14.0+2 

19.4+1 

q 7+1-0 
‘ -0.9 

10.6+1 

4.6 ±0.6 

8.1 ± 1.9 

488 

9 7+O.5 

Z- ‘ —0.5 

1090 

6-0+J 

13-4+1 

2 6 +0 ' 9 
z,o -0.8 

7-5+1 


Column 1: cluster name; Columns 2 & 3: projected mass within an aperture of 0.5//+ and 1/i+Mpc, resp.; Columns 4 & 6: va. (in units 


of h+kpc) determined using aperture masses; Columns 5 & 7: deprojected aperture masses within ; Columns 8-10: masses from best 


fit. NFW model All masses are listed in units of I I 0 171 1 1 VT. 1 











Updated CCCP weak tensing masses 33 


or the best fit NFW model. On the other hand, the estimate 
for JV /500 increases by 5% for the NFW model fit, whereas the 
deprojected aperture mass decreases by ~ 3% (also see §4.3 
in H12). Hence, to allow for a straightforward comparison 
with previous mass measurements, we present in Table ED 
the results i f we u se the mass-concentration relation from 
I Duffy et al.l J2008I '). 






